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ABSTRACT 


A conceptual design of a stand-off weapon to be launched 
from maritime patrol aircraft for use against hostile surface 
combatants was performed at the request of the Naval Air Test 
Center. The purpose of this thesis was to study the 
feasibility of developing a low-cost, anti-ship missile for 
alr ASW platforms. 

A mission threat analysis was conducted to determine the 
lethality of probable targets and to determine required 
missile performance characteristics. Current design methods 
and techniques were used to calculate the necessary missile 
geometry to meet the derived performance characteristics. 

An evaluation of navigation laws was conducted to 
determine the most appropriate flight profile for the missile. 
The control system was tallored to meet the specifications of 
the selected navigation law. 

An investigation of passive and active homing devices was 
conducted. A low cost seeker to adequately locate and track 
targets of interest was examined. 

A target engagement model was used to verify the missile's 
maneuverability. This model demonstrated that the missile 
could intercept highly maneuvering craft when launched from 


a desirable stand-off distance. 


Thala 


TABLE OF CONTENTS 


Mee INTRODUCTION 


A. 


ee 


D. HESTORICAL REV EERW gee meee. 6 oe eee 


E.« DESIGN GOADS Oo. i). roe cements meneen si.) ccememeenc caer 
Il. MISSION=THREAT ANADYS I See. 3s ee so eee 
A. INTRODUCTION a 0am ° So Coord Sicko et a6 oo 
B. DEVELOPMENT OF THE FRIGATE Class 
MODE 1s 2.53. 2 ae ° A eon oo Oe be 
C. DEVELOPMENT OF THE VSUBMARTNESCLASS 
MODEL... 5 202 ease eee ee So eet eo oe oe eer 
D. DEVELOPMENT OF THE PATROL CRAFT MODEL..... 


HISTORICAL REVIEW OF AIRPPROUECTE Ow. 


SUBMARINE WARFARE ¢i-02 2. o%. « ¢ eee ia aad 
MISSION DESCRIPTIChE. a < ateetheaenel. ee 
1. Forward Area ASW........ ee a 


OF Open Ocean ASW * ¢ @© @© @ @ oe @ @—6thmUhOHMmUC HhUCUc OhmhUCc OhUCUCcCmhUCc OhCUhhUh]}H HSH SF FS SF FF 


3. Direct ASW and SSSC Support of 


Naval FoOrceS..<<«e+2e2sees 6 pt eee 


4. Direct @SSBN/7SSN SUppoaE.-. ee. . «ee 


5. Direct ASW and SSSC Support of 


Merchant Shipman. ee 


6. Anti-Surface Warfare......02.)) ee 


7. Ocean Surveillance and Intelligence 


Opera EOns<ign cere BO a ar 
8. =Min Linge. . Veet eee o's nce ce (Pes eae ate 


JUSTIFPICAT LON. sesecne en eeecnchs cs ese ee 6 oem 


AY 


at Ee 


1a 


NAVIGATION LAW SELECTION........« SPORE es 6s 0 3 6 6 fu 
Ps TES) DIOS IR TOU 2 Sea 
B. CORRELATING THE NAVIGATION LAWS....... wens Wipe 
1. PROPORTIONAL NAVIGATION........ me secemer cts .< 
Care Omer NAVIGATION. 2... ct tte tet te ee 
Creeley ~DEARLNG NAVIGATION. .....-..ee.0-% 
e DETERMINATION OF REQUIRED LATERAL 
Pe GE iE eee Sos. a6. « 2 Wee of oe Sh 00 4 > ee 
eee iy Ors lg @ nWOIN chs « oss « scsQeMieMeisclieelWlsns « « clea se 0 « 
Pere ey ECOPMENT AND DES IGN 6 1 te 
io EN TRODUGRUON . «ices <u. fe, 
Boo HISTOR TGA LSSURY EY cumin. »-< S55 0 0 rr 
Gy ENITIAG SIZING OF THE MUNEGREONS CARRIER 
CLONE IG CRISS ELON Ss 6 i 
Popeorn ee TIONAL MISSILE ANALYSIS.......scecsecee 
Ee ee ESM Gle es % ¢ io see eo 6 0 «6 eisliee6@ 6's s+ 6 
1. Guidance/Seeker Weight............ - 
PeeCinemoat  BOdyY Siredl i. fi... .. 551 kt ee 
[Ne lhcwotesimgle Wing Panel......s.s...- 
F. DETERMINATION OF THE WING-BODY-CANARD 
ieee CURY BesGitODE cis... . 0 OM, ww ek 
G. DETERMINATION OF THE PITCHING MOMENT 
LOIS NS MIS EC) 2) SURE As ar arr On i 
H. HORIZONTAL TAIL-SIZING AS DETERMINED 


Sr oie ee DERECTIONAL STABILITY 
IS OME RSE) SES) 5 6 Ge Ceo. ca a mn ee 


Wy 


Vi 


Vat 


Li DRAG DETERMINATION Soe crt cememertemen. ss sc sentemercrs 
J. DETERMINATION OF THE PITCH AND YAW 
MOMENTS OF INERTIA< « « 2 5 sigier ence nen ne cee ee 
K. FINAL BASELINE CONE TCURAT 1@ Miscececon ene 
PROPULSION..... ee OO Ere ¢ 
A. INTRODUCTION tie oicus cette eee Sh “oi etents tae tate one 
Ea BOOSTER SLiZING were eee eee Se er ere hs 
C. SUSTAINER StZINGe. re Pe er 
D. NOZZLE SIZTNG ae PP aa es Sa on 
SEEKER DEVELOPMENT «.2e2..« 5 citemereitres) sete ct mee 
A. INTRODUCTION ee oa be See ueetenecke: 6 ohe cortece tet etremrerte ram 
B. APERTURE DETERMINATION .22...2) soe eteeene = cee nee 
1. Target Signature Estimation. ee 
2. IR Rangemhquation.....5. co 3 SR ee : 
C. SEEKER CONs LRUGTION sae rs Cnr ere ne 
By SEEKERWOPERAT LON « .0.ic) s:.cieweueneesiteieiemciteit Meiettctrc ron 
ES SEEKER EQUATIONSS OF MOTTON eee tl... soe 
AUTOPILOT DESIGN «2.42 ee Eo 
A. INTRODUCTION eee. ere of een A 
B. STABILITY. DERIVATIVES @oiviocuseeens « ate 
C. DEVELOPMENT OF THE AERODYNAMIC TRANSFER 
FUNCTIONS (ATF) = 0ec < c:c2ee: . e 
Ds LATERAL AUTOPILOT OPERATIONS types... se 
E. APPROXIMATIONS OF AUTOPILOT COMPONENTS.... 
F. AUTOPILOT STABILITY CONSITDERAWEG Noa... cm 


v0 


G. 


APPROXIMATION OF THE AUTOPILOT TRANSFER 


PWINGUUOChire ts sett ss ss se 6 es PROUT ss oe te wt 8 tem! 
Pie veeGH WME IIGAGEMENT OCENARIO.. 2... 1. ee te ee oo 154 
rx. Generis else D RECOMMENDATIONS. .. 2.2 cece er cece je Al 
Pe Ore MOE 6 6s ee se ee ss 8 es 8 Meee ss ee ss oo ss Por 
B. INITIALIZATION AND LAUNCH MODE 
DES GRal PIG Misa. swish. sis + 6s soso 6s s Ge. eye 4 6 ements 164 
Cech EKERS DESCRIPTION... ..2 Se 6 OO ee ee bey 
D ae COMMENDAT IONS... ....% Peal + siinkale co) ewememn. a cteme 166 
Bist OF REFERENCES..... SU eho ce ra hasa lo or sivaliaiiaapenwnaiis 169 
ne ree orem PON ELST a. fe ee ce ew ee NU 


Veal 


ACKNOWLEDGMENTS 


I am particularly indebted and grateful to my thesis 
advisor, Dr. R.M. Howard, for his professional counsel 
throughout the development of this thesis; his unwavering 
commitment to the success of this project was thoroughly 
appreciated. 

I would like to thank the following professors for their 
assistance in the evolution of this thesis. Thanks to Dr. 
H.A. Titus for his comprehensive review of this thesis and his 
notable suggestions in the development of the guidance system. 
Thanks to Dr. R.E. Ball for his expert advise concerningsene 
missile seeker design and the mission threat analysis. Thanks 
to Dr. E.M. Wu for his composite airframe structural analysis. 
Thanks to Dr. A.W. Cooper for his in depth explanation of 
infrared seeker operations. Thanks to Dr. A. Gerba, Dr. J. 
Hauser, and Dr. D.J. Collins for their suggestions on 
designing the autopilot system. Thanks to Dr. D.W. Netzer 
for his recommendations in the development of the solid rocket 
motor. And thanks to Dr. R. Wood for his insightful views of 
missile design. 

A special thanks to Mr. D. Knutsen of the Anti-Surface 
Planning Office for hosting my visit to NWC, China Lake. Mr. 
Knutsen's expertise provided essential information that made 


this project feasible. 


Vidi 


fewould Hrke €oO thank Ledr. C. Bess and Ledr. C. Holm for 
their assistance in developing the missile flight profile. 

Finally, I want to thank my most understanding wife, 
Nancy, for her affectionate support and love that allowed me 
to contribute my best effort to this thesis. Without her 
heartening encouragement and unselfish dedication, this thesis 
would not have been possible. To John Joseph, Victoria Marie, 
and Benjamin Michael, I thank you for your understanding and 


love that gave me the motivation to complete this work. 


1X 





l. INTRODUCTION 


A. HISTORICAL REVIEW OF AIR PROJECTED ANTI-SUBMARINE 
WARFARE 


During the early part of World War Two, when the German 
U-boat campaign was stripping the Allies of their much needed 
supplies, aircraft such as the Lockheed PBO-1 Hudson were 
entering service to neutralize this submarine scourge. The 
war effort also brought into service such aircraft as the U.S. 
Navy's Lockheed PV-1, the first aircraft designed specifically 
for anti-submarine warfare (ASW). The PV-1 extended convoy 
alr coverage from the east coast of the United States out to 
the mid-Atlantic sector which aided in diminishing the U- 
boats' effectiveness. Several types of aircraft were used for 
ASW during World War Two; those produced by Lockheed shall be 
examined here. This examination will show how maritime 
strategies had influenced Lockheed in designing the current 
Pemerol aircraft. (Ref. l:p. 2} 

The maritime strategy saetawee by the German High Command 
was to control the sea-lanes with a submarine force instead 
of a surface force. Until 1942 the U-boats were able to 
operate off the east coast of the United States essentially 


unscathed. Aircraft in the ASW role were extremely effective 


against the Axis submarine force. This competence was so 
devastating against the U-boats that the submarines were 
forced to operate in the mid-Atlantic region, safe from enemy 
air patrols. After this time, the U-boat campaign lost 
momentum and the battle for the Atlantic turned in favor of 
the Allies. [Ref. l:p. 2] 

The late 1950's saw little change in the role of patrol 
aviation. The primary mission of patrol remained ASW. ASW 
was generally practiced ina "neutral threat environment." 
(Ref. (2p) ce Cazenove defines the neutral threat 
environment as an environment in which the aircraft "would not 
draw hostile fire, although it was in the vicinity of a 
hostile platform, 1.e., a submarine." (Ref. 2:p. aie 
Therefore, little effort was put into the survivability 
aspects of the aircraft and the majority of the design work 
was devoted to improving the aircraft's reliability. 

During the same time frame, the Soviets were employing a 
"German Naval Policy" (Ref. 2:p. 22] which was to control the 
sea-lanes with a submarine force. Concurrently, the Lockheed 
P-3 Orion was in its developmental stages and was designed 
according to the perceived threat. Like past ASW aircraft, 
the P-3 was not expected to encounter hostile fire; the design 
reflected the goals of endurance and reliability. [Ref. 2: 


| Shel. 7] 


B. MISSION DESCRIPTION 


This section is based on the work of Cazenove [Ref. 2]. 
Unlike earlier patrol aircraft, the P-3 has acquired several 
mission roles outside ASW. Table 1 outlines the majority of 
the missions a P-3 crew would be expected to conduct. It can 
be seen from Table 1 that the P-3 has moved into a multi- 
mission role which exposes it to a greater number of hostile 
platforms. Each of these missions shall be described briefly 
and the threat environment defined. 

a Forward Area A8Ww 

Forward areas are defined as "those combat areas 
which are nearest to enemy concentrations and/or behind 
established enemy lines." (Ref. 2:p. 34] An ocean or sea 
located in a forward area is referred to as a transit lane. 
If hostilities were to erupt in Europe, an example of a 
transit lane would be the Norwegian Sea. The P-3's mission 
in these areas 1S to search, locate and destroy enemy 
submarines. 

The P-3's survivability is considered low due to the 
close proximity of operations near hostile territory. The 
"threat encounter probability (TEP)" [Ref. 2:p. 34] can be 
defined as the probability of an aircraft drawing hostile 
fire. In this case, the TEP is high because of the P-3's 


susceptibility to air and patrol boat attacks. 


Table 1. P-3C Operational Missions [Ref. 2:p. 33] 


Miss LeM 

Forward Area ASW 

Open Ocean ASW 

Direct ASW and Surface/Sub-surface 
Surveillance and Communication (SSSC) 
Support of Naval Forces 

Direct Nuclear Powered Submarine (SSN)/ 
Nuclear Powered Ballistic Missile 


Submarine (SSBN) Support 


Direct ASW and SSSC Support of Merchant 
Shipping 


Anti-Surface Warfare (ASUW) 


Ocean Surface Surveillance and 
Intelligence Operations 


Mining 


Threat Encounter Probability (TEP) Key: L-low 
M-medium 
H-high 


EH-extremely 
high 





Foe Open Ocean ASW 
The mission is the same as that for the forward area 
except that it is conducted in an open ocean scenario. The 


operating area is so vast that there is only a remote 


possibility of encountering hostile forces. Accordingly, the 
Hee is low. This mission profile is the type the P-3 was 
initially designed to conduct. 
3" Direct ASW and SSSC Support of Naval Forces 
This mission's requirements are to conduct ASW in 
cooperation with friendly naval forces, to supply surface 
surveillance as requested, and to provide communication 
support. The friendly units conduct operations in either the 
transit lanes or in an open ocean setting. Therefore, the TEP 
will vary with the area in which the mission is performed. 
4. Direct SSBN/SSN Support 
In this mission the P-3 will assist in localizing and 
destroying hostile targets that pose a threat to friendly 
submarines. The aircraft can, as well, provide communication 
relay when necessary. Due to the nature of an SSN's mission 
(ASW/ASUW), the P-3 will be faced with a medium to high TEP. 
5". Direct ASW and SSSC Support of Merchant Shipping 
The purpose of this escort mission is to increase the 
depth of a convoy's defensive zone. The mission requirements 
are the same as those for the Direct Support for Naval Forces. 
Because hostile forces will attempt strike operations against 
the convoy, the TEP is rated high. 
6. Anti-Surface Warfare 
Now equipped with the Harpoon anti-ship missile 


system, the P-3 has become a formidable adversary to enemy 


platforms. The TEP is dependent on whether the crew can 
employ the Harpoon's over-the-horizon capability. If this 
capability can be utilized, then the TEP is rated low. If, 
however, the target must be identified visually prior to 


firing, the aircraft's survivability will be drastically 


reduced. The TEP would be rated high for the latter 
Situation. 
Th Ocean Surveillance and Intelligence Operations 


This mission requires the aircraft to search an area 
of ocean, identify, and track assigned targets. If the target 
is of special interest, intelligence information can be 
collected with sensors on board the P-3. 

If this mission is conducted prior to hostilities the 
TEP can be considered low. However, conducting such a mission 
in a hostile environment would increase the TEP significantly. 

8. Mining 


As the name implies, the P-3 would be required to 


drop mines in forward areas of interest, i.e., a harbor or 
transit lane. This mission normally requires more than one 
aircraft. A formation has a higher probability of being 


detected than a single aircraft. Consequently, the enemy may 
have an opportunity to repel the mining intrusion. The TEP 
for this scenario is rated high and approaches unity if the 
enemy is alerted. This mission is considered extremely 


hazardous for the P-3. 


Because the P-3 is a highly vulnerable aircraft, even 
a low TEP rating may prove fatal. Reviewing the missions and 
taking into consideration the TEP for each, one can conclude 
that the P-3 would have a potentially low survivability rate. 

A low survivability rate can imply an excessive loss 
rate. There are two options available to counter an excessive 
loss rate. The platform can be produced in such immense 
numbers that a high loss rate can be tolerated, or the 
platform can be modified to enhance its’ survivability 
characteristics. Mass production may not be possible owing 
to construction time and cost per platform, as could be the 
case with the P-3 Orion. Therefore improving survivability 
would prove to be a more viable alternative. Improving the 
survivability rate would require several courses of action. 
First, the P-3 must be eliminated from missions that would be 
considered extremely hazardous. This action would be the 
least expensive, but would unfavorably reduce the 
effectiveness of the platform and limit its areas of 
operation. Modifying the aircraft's defensive abilities may 
provide a compatible solution for survivability enhancement. 
Possible modifications are airframe improvements to reduce the 
aircraft's vulnerability, the addition of sensors that would 
alert the crew to an impending attack, and an expanded weapons 
inventory. The addition to the weapons inventory should 


include air to air missiles and air to surface missiles (ASM) 


that can be used in a defensive manner allowing the crew to 


disrupt the attacker's firing solution... (Ref. 2:p. 49) 


C. JUSTIFICATION 


Naval Air Systems Command (NAVAIRSYSCOM) has tasked the 
Naval Air Test Center (NATC) to investigate feasible solutions 
for improving the P-3s' survivability rate. One of these 
solutions is the addition of a short range ASM to the P-3 
weapon inventory. Table 2 specifies the desired parameters 
of such a weapon. According to NATC, current inventory weapon 
systems are unsuitable for the patrol (VP) mission profile. 
The Naval Postgraduate School has been requested by NATC to 
assist in the development of a stand-off weapon for patrol 
aviation use. (Ref. 3] 


The purpose of this study is to conceptually design a 


weapon to meet the parameters outlined by NATC. The Thesis 
is divided inte five parts. Part one develops required 
background information. In Part two a threat analysis is 


conducted. From this analysis, target characteristics can be 
developed into models for simulation use. Part three develops 
a guidance plan based on the models produced earlier. Part 
four designs the propulsion system, airframe, seeker head, 


and control system. Part five is a cost analysis and summary. 


Table 2. Stand-Off Weapon Design Parameters [Ref. 3] 


RANGE: 0-15 NM when launched from 1000 ft. AGL 


PERFORMANCE*#: Launchable from P-3 service ceiling; 
allow aircraft to begin evasive 
maneuvering 












GUIDANCE: "Fire and Forget" capability 


Patrol craft 





PERSPECTIVE TARGET SIZE: 





Frigates (under 1K TON) 






Surfaced diesel submarines 






WARM-UP TIME: Short as possible so weapon may be used 
in a defensive nature as well as 
offensive 


COST: Less than $200,000 per weapon 
* Weapon should be compatible with current armament 
system. 


D. HISTORICAL REVIEW 











A consideration of selected ASM used by patrol aircraft 
will assist in avoiding the imperfections experienced in past 
designs. This section was developed from Reference 4. 

The origins of the ASM began in 1937 when RCA became 
interested in a television guided weapon system. During 
August 1940, RCA proposed a television guided glide bomb to 
the National Defense Research Committee (NDRC). By January 
1941, a glider with a 12-foot wing span capable of carrying 


a 2000-pound bomb had been developed, designated "Dryden". 


In April 1942, the Army had developed a semi-active radar 
guided bomb. The semi-active seeker head was matched with the 
Dryden airframe to produce the glide bomb designated 
"Pelican", which was fitted on the Navy's PV-1 aircraft. 
Unfortunately, the operators had difficulty locking on to 
targets and the operational range of the PV-1 had been reduced 
by 20% with the addition of the glide weapon. Because of 
these problems, Admiral King cancelled the program in 
September 1944. During the same time period, an active radar 
guided glide bomb was under development. This weapon, 
designated "Bat", used the Dryden airframe carrying a reduced 
payload of one 500-pound general purpose bomb. Bombing 
through the overcast (BTO) radar brought the necessary 
technology needed to produce the Bat's seeker head. The Bat 
was the first ASM to enter combat. By May 1945, the Bat had 
scored three direct hits against enemy ships, and had nearly 
a 40% hit record. The advantage of the Bat missile is that 
it allowed the attacking aircraft to remain outside the lethal 
range of the anti-aircraft batteries aboard hostile vessels. 

The delivery method involved the launching aircraft 
acquiring the target on the missile's seeker head, then 
releasing the weapon at a range of 15 to 20 nautical miles. 
There was a drawback to these comfortable ranges; the aircraft 


needed a launch altitude of 25,000 feet, which left the 
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Gelivery aircraft open to fighter attack. Therefore, Bat 
delivery aircraft usually had fighter escort for protection. 

During post-war testing, the Bat was found to have a 
serious flaw; the radar seeker would become saturated if other 
radars were in operation near by. Unable to resolve the 
problem, the program was cancelled in 1948. 

When the P-3 Orion entered service in the early 1960's, 
the Soviets had armed designated classes of submarines with 
long range surface to surface missiles. These missiles, 
designated SSN-3 by NATO, were to be utilized in anti-carrier 
operations. The P-3 was equipped with the air to surface 
Bullpup missile as a deterrent to the missile equipped 
submarine menace. 

The Bullpup was a radio-controlled supersonic ASM. 
Carrying a 250 pound warhead, the weapon could be launched at 
a maximum range of Seven miles. The steering commands were 
generated by the P-3 co-pilot, lining the target up through 
a gun sight; he would make course corrections by observing a 
flare mounted in the missile's tail. 

A fault of the Bullpup ASM was its limited range. The 
P-3 aircraft had the capability to espy a surfaced submarine 
at long distances. In spite of this capability, by the time 
the aircraft got within range to fire the Bullpup, the 


submarine could have launched its missiles and submerged. 
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In 1970 NAVAIRSYSCOM released a Request For Proposal (RFP) 


for an ASM that would give patrol aircraft a long range strike 


capability against surfaced submarines. Mc Donnell Douglas 
was selected as prime contractor in 1971. The name given to 
this ASM was "Harpoon". The name was derived from the 
mission; "a 'harpoon' to attack 'whales.'" [{Ref. 4:p. 211] 


The Harpoon was originally designed to destroy the 
pressure hulls of submarines with its 510-pound penetrating 
blast warhead. Harpoon was soon adopted for ASUW as well. 

The guidance system is of the "fire and forget" scheme. 
Once programmed with targeting information, the Harpoon 
requires no further link with the launching platform. Ata 
predetermined range, the frequency agile radar initiates a 
selectable search pattern for the designated target. Upon 
target detection, the Harpoon drops to a sea-skimming profile. 
When the missile approaches the terminal phase, it performs 
a pop-up maneuver and dives into the target. 

"Ranges as great as 60 nm have been reported" [Ref. 4:p. 
212), which give the P-3 a significant stand-off capacity. 
The extended range capability also requireS a more 
sophisticated weapon system which drives up the weapon's cost. 
Each Harpoon is a 1.6 million dollar expenditure [{Ref. 5: p. 


183]. Cost is Harpoon's major drawback. 


de 


E. DESIGN GOALS 


Before the actual designing commenced, primary design 
goals were established. In the case of the stand-off weapon 
venture, NATC has provided a set of parameters which are to 
a large extent flexible. 

The weapon must possess the "fire and forget" scheme. 
This scheme allows the crew to position the aircraft ina less 
susceptible locality. 

Designing a weapon to give the launching platform a 
defensive capability appears uncommon to ASM. Usually, ASMs 
are utilized in an offensive manner, implying the attacker 
attempts to approach the target undetected. A typical 
offensive strike may proceed as follows: the target is 
located, targeting information is fed to the missile, and the 
incursion is enacted. In a defensive role, however, the 
Meeerle platform is already in a firing position, and the 
friendly unit must maneuver immediately to avoid being hit. 
If certain weapons onboard the friendly unit did not require 
targeting information, the crew may have the opportunity to 
deliver an immediate Pountersatracn against the firing 
pratiLorm. 

As an illustration of the above conditions consider a 
Simplified version of : Maritime patrol ASM strike. Assume 


a target is seen by the pilot, the target's position is 
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evaluated and placed in a format that can be utilized as 
targeting data. These data are inputted to the ASM and the 
weapon is launched. In an offensive situation this routine 
is both adequate and necessary to insure the highest 
probability of a kill. However, the extra time spent by the 
crew to program the missile in a defensive scenario may prove 
to be catastrophic. 

If a weapon were accessible with a "point and shoot" 
capability, missile programming would not be required. The 
pilot, upon sighting the target, could point the air¢Grameeeen 
the direction of the target and release the weapon. This 
capability would give the crew the ability to deliver an 
immediate counter-attack, perhaps deterring any further 
attacks from the enemy vessel. 

The goals listed above are driven by cost; an elaborate 
system would probably exceed the allowable price per weapon 
Maguire: Striving to keep costs down, "off-the-shelf" 
components will be used where feasible. Using this procedure 
removes research and development expenditures for individual 
components. Another cost reduction scheme would incorporate 
utilizing the existing armament delivery system presently 
onboard the P-3. Employing these techniques will facilitate 


in preserving the goal of a low cost weapon. 
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ll. MISSION-THREAT ANALYSIS 


A. INTRODUCTION 


The purpose of this analysis is to determine the required 
missile performance parameters necessary to impose a kill 
against specified targets. This analysis will consider the 
targets' maneuverabilities and anti-air capabilities. The 
threat posed by the target against the delivery aircraft is 
not considered in this chapter. The missile system under 
development is designed to keep the aircraft outside the 
lethal range of anti-air weaponry on board target platforms. 

As required by NATC, the target size must range from a 
small frigate (less than 1 kiloton) to a surfaced submarine. 
A variety of platforms are encompassed in this size range. 
Therefore, target models will be developed to generalize the 
threat characteristics of those platforms incorporated in the 
target size parameter. These models will be selected from 
three classes of ships. The designated classes are submarine, 
frigate, and patrol craft. Each class represents a unique set 
of characteristics which needs to be considered in the design 
of the stand-off missile. The submarine class depicts a slow 
moving platform when surfaced, with a structural integrity 


that resists missile penetration. The patrol craft displays 


Its, 


the ability to maneuver rapidly, which forces the missile to 
develop extreme lateral accelerations to achieve an intercept. 
The frigate exhibits moderate maneuverability, but the primary 
threat to an incoming missile is the frigate's anti-air 
defense systems. 

An analysis of existing platforms in each of the three 
model classes waS conducted to determine target threat 
characteristics. These characteristics were used in the 
development of the target models. The approach used to 
conduct this analysis was to examine several navies of the 
Warsaw Pact and navies of those nations receiving military aid 
fromthe Sov Leteuniion. A nation was considered for this 
analysis if it had received over 100 million dollars in Soviet 
arms over the period 1981-1986. Table 3 lists the countries 
meeting this criterion and the dollar value of Soviet arms 
transferred. The navies of these countries were surveyed for 
the most common platform in each model class. The 
characteristics of these common platforms would constitute the 
parameters needed to develop the generalized models. Table 
4 recapitulates the countries examined, the total number of 
platforms in a given model category, and the class name of the 


most prominent platform for each model category. 
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Table 3. Exported Soviet Arms in Dollars [Ref. 6:p. 22] 


DOLLARS RECEIVED COUNTRY 


57 LO LOsSLELION 


1 TO 5 BILLION ALGERIA 
IRAN 

100 MILLION TO 

1 BILLION NICARAGUA 


B. DEVELOPMENT OF THE FRIGATE CLASS MODEL 











The predominant platform in the frigate class is the 
Soviet Grisha III. The Grisha III is a small anti-submarine 
ship ("Soviet type designation: MALYY PROTIVOLODOCHNYY KORABL" 
(Ref. 7:p. 535]). The parameters of the Grisha III are listed 
in Table 5. The armament system of the Grisha III provide it 
With an impressive anti-air capability. This attribute has 
the potential to reduce missile survivability; therefore, the 
Grisha III armament system shall be incorporated into the 
general frigate model. 

The frigate class model developed from the analysis has 
its parameters listed in Table 6. Figure 1 shows a line 
drawing of the model and Figure 2 depicts the lethal ranges 


of the ship's anti-air weapon systems. 


ily 


Table 4. Prominent Platform for Model Category [Ref. 8] 


COUNTRY FRIGATES PATROL SUBMARINES 
CRAFT 


SOV LET rege TOT 275-161 
UNION 17% GRISHA III 29% OSA 


CUBA 2 TOT _60 TOT 
100% KONI 30% OSA 


LIBYA 2 1enk 
100% KONI 







147 TOT 
40% FOXTROT 











SLOT 
100% FOXTROT 






6 Ter 
100%. FOxXTRGs 











13 er 
COMBATTANTE 











SYRIA 22 LOT Zee Low 
100% PETYA 44% OSA 
ALGERIA 3. LOm 22. TOM 2-1ot 









100% ROMEO 


100% KONI 92% OSA 





TOTAL NUMBER OF SHIPS IN CATEGORY 
PERCENTAGE OF TOTAL/CLASS NAME 
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Table 5. Platform Parameters 


| GRISHA III | FOXTROT | OSA I! | 
(Ref. 8:p. 580) (Ref. 7:p. 580} (Ref. 7:p. 542) 








| 
| DISPLACEMENT 950 TONS 1950 TONS 215 TONS 
MAX SPEED 15.5KTS (SURFACED) 


DIMENSIONS (HULL) | 236 X 33 X 12 FT 297 X 23 X 20 FT 126 Xe2aeXe 7 ET 





ANTI-AIR ARMAMENT | 2-SAN-4 LAUNCHERS 2-TWIN 30MM 
GUNMOUNTS 
1-30MM GATLING | 
| MACHINE GUN | 
| 1-TWIN 57MM | | 
| | | 


| GUNMOUNT 
ce 


C. DEVELOPMENT OF THE SUBMARINE CLASS MODEL 


Examinations of Table 4 reveals the Foxtrot, a Soviet 
diesel submarine, to be the most prevalent platform in this 
category. The armament system of the Foxtrot appears to lack 
an anti-air capability; however, the Soviet exported Kilo 
class diesel submarine seems to possess some type of SAM 
system. (Rei. 7.p. 507) #&Thus’, the technology appears to be 
available for a submarine to protect itself from an air 
threat. Table 5 lists the Foxtrot's parameters. The 
submarine's limited ability to maneuver on the surface is 
assumed to pose little challenge to the control system of an 
air-to-surface missile. The submarine, nevertheless, presents 


two distinct quandaries. First, the submarine has the ability 
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Table 6. Target Model Parameters [Ref. 9] 


| FRIGATE SUBMARINE PATROL CRAFT | 


DISPLACEMENT 950 TONS 2000 TONS 215 TONS 
MAXIMUM SPEED 30KTS 16KTS (SURFACED) 35KTS | 


DIMENSIONS 235 X 33 X 12 FT 3°70 X 26 X 20 FT 130 X 26 X 7 FT 
| 
ARMAMENT 2-SAM LAUNCHERS 1-SAM LAUNCHER 2-TWIN 30MM MOUNTS 
RNG: 8NM RNG: 1. 75NM RNG: 1.6NM 
PH: 49% PH: 18% PH:21% 
1-30MM GATLING GUN 
RNG: 1.5NM 
PH: 42% 
| 
| 1-TWIN S7MM GUN | 
RNG:2.7NM 
PH:35% 
| MAX. HEADING 
CHANGE leo°7s 0.33°/s 6°/s 


TURN RADIUS 0. SNM i 0.250 
MAX. ACCEL. 0.05KTS/s 0.03KTS/s 0.16KTS/s 


| KEY: PH - PROBABILITY OF HIT RNG - RANGE MAX. ACCEL. - MAXIMUM ACCELERATION | 
a a 


to submerge. If the submarine's crew became aware of an 
impending surface attack, they could remove themselves as a 
surface target. The second problem is the submarine's double 
pressure hull construction. Reference 8 cites that "the 
patrol class submarine is similar in design to the German U- 
boat type XXI." [{Ref. 8:p. 551] The type XXI had a deuamte 
hull construction and the hulls "were formed carbon steel 
plating 28mm thick." [Ref. 10:p. 76] The thickness of these 
pressure hulls may reduce the effectiveness of the missile 


warhead. Possibly, the warhead may be able to create shock 
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see Table 6 for dimensions 


~— 
a) ; 
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Frigate 





| | 
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Submarine 
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Patrot Cratt 





Figure 1. Target Model Profiles. 
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see Table 6 for range data 
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| 
Frigate 

| 

| 

| Submarine 
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Figure 2. Target Model Range Rings. [Ref. 9] 
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damage to the submarine's structure so it would force the 
submarine to remain surfaced. Such shock damage brought a 
Soviet Yankee class SSBN to the surface off the coast of 
Bermuda in 1986. Missile fuel ignited, causing damage to the 
Yankee's structure [Ref. 1ll:p. 12}. If a submarine is unable 
to submerge, its usefulness as an offensive weapon has been 
greatly reduced. 

Table 6 displays the parameters for the submarine class 


model. Figure 1 shows a line drawing of the model. 


D. DEVELOPMENT OF THE PATROL CRAFT MODEL 


The most common patrol craft in the analysis is the Soviet 
Osa class. The Soviet designation for this class is "RAKETNYY 
KATER". {Ref 7:p. 542] The Osa's main offensive weapon, the 
surface-to-surface STYX missile, can pose a serious threat to 
surface units. These patrol craft also possess an air defense 
capability, which could hinder an air strike using free fall 
weapons. The weaponry and other specifications of the Osa 
are displayed in Table 5. Some units are reported to be 
fitted with the surface-to-air missile system, which is a IR- 
homing, visually aimed anti-air system. [Ref. 7:p. 542] 

Besides the craft's anti-air systems, it possesses another 
obstacle for the guidance system of an incoming missile. 


These small craft are highly maneuverable; thus the guidance 
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system must continually adjust the missile's position to 
achieve a successful intercept. 

Figure 1 shows a line drawing of the patrol craft model 
and its parameters are displayed in Table 6. 

The models produced in this chapter will be used to select 
the proper navigation law, to determine warhead lethality, and 
to estimate the probability of a kill against a designated 


Carder. 


24 


Il. NAVIGATION LAW SELECTION 


A. INTRODUCTION 


A navigation law is defined as "the analytical formulation 
used by the guidance system to convert sensed target 
information into missile steering commands." [Ref. 12:p. 35] 
The objective of a navigation law is to determine the 
necessary lateral acceleration the missile need generate to 
achieve an intercept with the target. Lateral acceleration 
is developed by the lifting surfaces of the missile. 
Consequently, the airframe is designed to satisfy the lift 
condition needed to produce the perceived required lateral 
acceleration. 

There are four general navigation laws. They are: 

fo Pursuit 

- Constant Bearing 

+ Line-of-Sight (LOS) 

- Proportional Navigation 
Pursuit and constant bearing are special cases of proportional 
navigation. Their similarities will be discussed later. LOS 
is used primarily for beam riding weapons. These weapons 
normally require a tracker to be mounted on the delivery 


platform. The tracker illuminates the target so the missile 
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will home on the reflected energy. LOS, however, does not 
meet the required "fire and forget" capability because the 
launching platform must track the target for the missile. 
Hence, only proportional navigation and its derivatives will 
be explored. (Ref. 132). 3674) 

Pursuit navigation is "a course in which the missile 
velocity vector is always directed toward the instantaneous 
target position." [Ref. 14:p. 460] The advantage of pursuit 
navigation is that the navigation information is simple, which 
makes the avionics light weight and less expensive than for 
other navigation laws. Another advantage is, the missile does 
little maneuvering until it is close to the target's position. 
This advantage decreases the induced drag the missile would 
produce during the cruise portion of an engagement. The main 
drawback pursuit navigation has 1s encountered during the 
terminal portion of the flight: the weapon is in such a 
position relative to the target that the missile requires 
large lateral accelerations to make the target intercept. 
Pursuit navigation seems best employed against stationary 
targets. LRG fies ome] 

Constant bearing navigation is a "course in which the LOS 
from the missile to the target maintains a constant direction 
in space." [Ref. 14:p. 473] If the target has a constant 
velocity then the missile's required lateral acceleration will 


be zero. If the target maneuvers, then the missile's required 
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lateral acceleration will never exceed the target lateral 
acceleration. This advantage gives the missile the ability 
to intercept the target with minimal required lateral 
acceleration. However, the guidance must possess the ability 
to predict the future target positions. This necessity 
complicates navigation avionics, increases missile weight, and 
elevates system cost. 

Proportional navigation (pronav) is a "course in which the 
rate of change of missile heading is directly proportional to 
the rate of rotation of the LOS from the missile to the 
target" [(Ref. 14:p. 475]. The advantage of pronav is the 
increased sensitivity of the navigation system to target 
mobility. Pronav positions the missile in such a manner that 
during the terminal phase the required turn rate will be near 
zero. Pronav has the most complex avionics system of the four 
general navigation laws. This complexity implies extra weight 
and increased cost affixed to the missile design. [Ref. 13:p. 


67) 


B. CORRELATING THE NAVIGATION LAWS 


This section was developed from Reference 13; its 
objective is to show that correlating the navigation laws will 
alleviate the need to develop separate computer programs for 
each law. The geometries of the previously discussed laws 


are depicted in Figure 3. 


Z7 





- én a 
| y ZA a \ 


Constant Bearing 





PRONAV 


Figure 3. Geometries of the Navigation Laws [Ref. 13] 
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The above navigation laws can be universalized by the 
following statement: 


The missile rate of turn is a multiple of the rate of turn 
of the LOS, or in equation forn, 


@= KD one 


The multiple is known as the navigation constant, K. 
meet. 13:p. 67) 


Inspection of Figure 3 produces the following angle 


relationships for each law addressed: 


* PURSUIT p = ee (2) 

- CONSTANT BEARING He O ee ce (3) 

- PROPORTIONAL oP D - aa (4) 
Rearranging equation 1 to solve for Z Peer 


2= 6, / K (5) 
Equations 2, 3, 4, and 5 may now be used to_- show the 


correlation between navigation laws. 


i. PROPORTIONAL NAVIGATION 
OSES EE (6) 
Therefore, 
e L € 
= — q 
On= L- Ky (7) 
Miescituting equation 5 into equation 7: 
: \-K \ | 
eo (8) 
Oc. | Gyn 


2 PURSUIT NAVIGATION 
For pursuit navigation K=1, and substituting this 


value into equation 8 yields: 


Zo 


O m= O 
Therefore equation 6 becomes: 
Don= Gry 
which is equivalent to equation 2. 
a CONSTANT BEARING NAVIGATION 
e 
The definition of constant bearing implies D = OQ, 


provided the target maintains a constant velocity. Hence, 


equation 7 bears ; 


G. Of 
upon substitution of this value. This agrees with equation 


3 and suggests K must be selected to meet the equality of 


far) * | 


This section demonstrates that a navigation law is dependent 


equation 8. Therefore, 


on the value of K adopted. This similarity is useful in the 
development of a computer program that can be used to 
determine the foremost navigation law for the given target 
intercept scenario. [Ret ise pp. J6Go) 
There are other parameters which also effect the 

target intercept. They are: 

* Target Heading Changes 

- Target Speed 

-* Target Acceleration 


* Seeker Angle Bias 
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+ Seeker Noise 
+ Wind Gust Effects 
The effects of these parameters on target intercept 
have been analyzed by Goodstein {Ref. 15]. Presented in the 
analysis are several graphs and tables that provide an 
indication of the navigation laws' sensitivity to a given 
parameter. Constant bearing navigation was not examined in 
Goodstein's analysis. Regardless, the information that was 
Peesented on the remaining laws is sufficient for the purpose 
of a preliminary navigation law selection. Table 7 displays 
a simplified version of the study. The table shows how 
adequately a particular navigation law could compensate for 
a given parameter. After examination of Table 7, it is clear 
that pronav has superior qualities when compared to pursuit 
navigation; however, Goodstein [Ref. 15] comments that prior 
to the final navigation law selection, cost and complexities 
of the required avionics must be considered. (Ref. 12:p. 38) 


Table 7. Navigation Law Sensitivity Parameters [Ref. 15] 


| | TARGET | TARGET | TARGET | SEEKER | SEEKER | WIND | 


HDG CHG SPEED | ACCEL | ANGLE BIAS; NOISE | GUSTS 


l 4 


PURSUIT | PooR | = G000 | GOO0 








| PROPORTIONAL | GOOO 


ot 


C. DETERMINATION OF REQUIRED LATERAL ACCELERATIONS 


Modification of a BASIC program developed by Redmon [Ref. 
12:p. 68} is used to determine the required lateral 
acceleration needed to intercept the model target. 

The patrol craft class was the selected target model. 
This class exhibits the greatest maneuverability which poses 
the most challenge for a missile guidance system. 

In this analysis the lateral or normal acceleration 
produced by the target is of concern. The target's speed is 
considered constant throughout the scenario; thus, upon 
maneuvering, the tangential acceleration component is zero. 
The normal acceleration of the target during a maneuver may 


be calculated from: 
= oy 
hi a ee, (9) 
If 


where, 
ris the radius of turn 
vV is the speed of the target 
The patrol craft has a turn radius of 0.25 miles and a 
speed of 35kts as can be drawn from Table 6. It is assumed 
the missile is launched at a crossing target. The maneuvering 
target's lateral acceleration may be developed by substituting 


the patrol craft's parameters into equation 9. This produces 


BZ 


a target acceleration of: 
Q 


VA 


vs? 
(18)°/ 45.8 
Outwe/.oe 


N 


{! 


Table 8 presents the lateral acceleration determined by 
the computer simulation. Pursuit navigation, even under these 
simple conditions, demands substantial lateral accelerations 
during the terminal phase of flight. 


Table 8. Lateral Accelerations Developed by Navigation Laws 
BOR Ss Ui CONSTANT BEARING PRONAV 
ee 1) ae Se 10) 2 “oe 
LATERAL 
ACCELERATION 269 /S 0.6 m/s* 0.9 ny S 4 
DEVELOPED 


Figures 4, 5, and 6 show the path taken by the target and 











missile for each of the navigation laws discussed at a range 
of three miles (5.5 km). 

The general conclusion of this analysis is that pursuit 
navigation does the majority of its maneuvering during the 
last few seconds of an encounter, whereas constant bearing and 
pronav produce large maneuvers initially so that only small 
corrections to the missile's flight path are needed near the 


Semeclusion of the encounter. 
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Pronav Trace 
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Figure 6. Pursuit Nav Trace 
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The navigation law selected for this design is pronav. 
The qualities illustrated throughout this chapter make it the 
optimal choice for a maneuvering target scenario. The reduced 
lateral accelerations required lessen the demands on the 
airframe, which implies a lighter structure. The increased 
avionics weight of a pronav system would most likely be less 
compared to the increase in airframe weight if pursuit 
navigation was selected. If the pronav system is more likely 
to intercept the target than a system employing pursuit 
navigation, then the number of missiles required to sink a 
desired target could conceivably be less than that required 
by the pursuit navigation system. Hence, the cost per kill 
for the pronav system would be substantially lower than that 


of the pursuit navigation system. 
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IV. WARHEAD SELECTION 


If a powered munitions carrier configuration were selected 
for the overall missile design, the missile could be fitted 
with a payload that best suits the mission requirements. A 
design feature such as this enhances missile operational 
flexibility. 

The munitions carrier will be restricted to the 1000-pound 
weight class. This limit permits the missile to complement 
other weapon loads carried by maritime patrol aircraft. 
Therefore, the payload that can be outfitted on the carrier 
must be limited to 500 pounds. For the purpose of this 
design, the baseline payload will be a generic 500-pound 
general purpose bomb. 

The 500-pound bomb is shown in Figure 7. The bomb has a 
cast steel case and is loaded with 192 pounds of H6 high 
explosives [Ref. 16:p. 247]. Figure 7 also shows’ bomb 
dimensions and center of gravity location. 

If the conventional missile configuration is selected, the 
damage mechanism would be a shaped charge warhead. Two charge 
Sizes will be examined. These sizes are: 

¢* 150-pound warhead 


* 300-pound warhead 
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Figure 7. General Purpose 500-Pound Bomb [Ref. 17] 





If the missile can be tailored to fit in the 500-pound weight 
class, then it may be carried on any P=-3 wing station. This 
advantage allows the missile to complement other weapon 
configurations. Therefore, the warhead size will depend upon 
the weight of other components so the missile gross weight 
does not exceed the weight limit on any one wing station. 

A cutaway of a generic shaped charge is shown in Figure 
8. As can be seen, approximately 40% of the warhead is HE, 
while the remainder is dedicated to fuzing and structural 
support [Ref. 18:p. 19]. Lindsey [Ref. 13:p. 190] describes 
the warhead detonation sequence: 

When the shaped charge warhead strikes a target, the 


point-detonating nose-fuze fires a length of detonating 
cord which leads to a booster in the rear of the warhead. 


a9 


The booster in turn detonates the main charge and a 
detonation wave travels forward causing the metal front 
liner to collapse. Collapse of the liner starts at the 
apex. When the liner collapses, it ejects a narrow jet 
of explosive products and metal particles from the face 
of the liner out the front end of the thick casing at 
velocities from 10,000 to 38,000 feet per second. 


This series of events produces severe internal damage to the 


target. 


detonating chord nose-fuze 


liner apex 





Figure 8. Shaped Charge Cut Away [Ref. 13] 

The conversions introduced by Bond [Ref. 9:p. 56] (agi 
used to evaluate the effectiveness of the generic weapons when 
employed against the target models. 

The amount of damage a target model is capable of 
withstanding is related to its displacement. The displacement 


of each class will be converted into a point value. This 
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Table 9. Displacement Points [Ref. 9] 


—_Taisccnenie sonia 


point value is equivalent to the damage required to sink the 










target. Table 9 displays the point rating for each model 
class. 

For evaluation purposes, the destructive force of the 
generic weapons should be scaled with the same point system 
as the target models. 

The amount of damage each weapon can deliver is shown in 
Table 10. To calculate the required number of weapons needed 
to sink a particular target, divide the target's displacement 
Value by the destructive point value of the bomb. The minimum 
number of weapon hits required to sink the target is displayed 
Mmmerable 11. ([{Ref. 9:p. 56] 

In summary, the 500-pound and 150-pound shaped charge 
appear to produce the same damage effects; they seem to 
require numerous target hits to produce a kill. If the kill 
level is relaxed, where the damage incurred by the target 
would force it to head for the nearest friendly port; then, 
perhaps the 150-pound charge would prove to be the best 


selection for a conventional missile design. This warhead 
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Table 10. Single Weapon Damage Points [Ref. 9] 


500-POUND 

GENERAL PURPOSE l/, 
150-POUND 
SHAPED CHARGE ilies 
300-POUND 
SHAPED CHARGE PG | 


Table 11. Minimum Number of Weapons Required to Sink Model 
fRef. 9] 


| Friars | FRIGATE SUBMARINE | PATROL Ce CRAFT 
500-POUND 
GENERAL PURPOSE 
Va0>Peun 
SHAPED CHARGE 3 3 dL 


S00=POUND 
SHAPED GHA Ree 





















size would keep the weapon in the 500-pound weight category. 
The 300-pound shaped charge may prove optimal if instead of 
the kill level being reduced, the target engagement range were 
decreased. This reduced range may afford a lighter motor so 
the extra warhead weight could be accommodated; but from a 
historical review, most missiles with 300-pound warheads 
usually weigh close to 1000 pounds. Thus, for the 


conventional design the 150-pound shaped charge will be 
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Pemected and, for the munitions carrier, a 500-pound general 


purpose bomb appears to produce an acceptable kill level. 
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V. AIRFRAME DEVELOPMENT AND DESIGN 


A. INTRODUCTION 


Methods used to select and design the wing, body, and tail 
surfaces will be addressed in this chapter. A reader 
interested in the final airframe design is referred to Chapter 
9, Conclusions and Recommendations. 

There are two configurations under examination. The first 
configuration is classified as a munitions carrier. This 
design allows the missile external payload to be changed as 
dictated by the operational situation. The second design has 
the profile of a conventional missile with a non removable 
payload. The munitions carrier requires the accommodation of 
a five hundred pound payload. Chapter 4, Warhead Selection, 
discusses the payloads that are available for this selection. 

The conventional missile configuration requires the 
ability to carry a maximum payload of one hundred and fifty 
pounds. 

Both configurations are limited to a wing span of 34 
inches. This restriction is based on the minimum distance 
between wing weapon stations as mounted on the P-3 Orion. To 
keep the missile design uncomplicated, the used of a wing fold 


mechanism will be avoided. Accordingly, the selected design 
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must be able to sustain flight with a fixed wing span of 34 


miches. 


B. HISTORICAL SURVEY 


A historical survey was conducted to assist in the initial 
sizing of the missile components. Table 12 lists the 
pertinent information needed to size the configurations under 
consideration. 

Table 12. Missile Historical Parameters [Refs. 4,5,7,8] 


wee 





AR. | AR. | L/b [Mo/hn | tw/? aes [t/t Vol 





| MISSILE | L/D W/S |S /S ray?) 


| RB OSE | 8.95 | 68.3 [12.8] 2.4]7.5 | 2.3] - — | a 


GBv15 | 8.6 | 96.6 | 7.9 | 0.8 | 1.3 | 2.7] 2.6 | - : : 


RBS 15F 055 113.S0lmie2 | 7.5 \e3.1 | 





SEA SKUA 
PENGUIN GAS 0 \ 1 29UORS| 243 | 322° | 2.9 
MAVERICK 
SKY SHARK | area 


| TOMAHAWK | 12.1 | 164.4 | - | 6.7 | 















Subscript Key: M:Motor WH: Warhead AV: Avionics C:Canard 
| W:Wing WH Vol:Warhead Volume 


C. INITIAL SIZING OF THE MUNITIONS CARRIER CONFIGURATION 


This section will develop a baseline configuration from 
which a more detailed design effort may be initiated. The 
procedure used was developed by Redmon [Ref. 12:p. 143}. 


Table 12 is referenced throughout this section. 
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A canard configuration has been selected for this design. 
This selection allows major components to be placed in 
separate compartments. For example, fin actuators are not 
located in the aft portion of the motor section as is the case 
with an aft controlled missile. Such a configuration would 
require some form of interface to be routed through the 
warhead and motor sections to the actuators. This interface 
complicates missile construction and maintenance efforts as 
well as drives up cost. [Ref. 19:p. 45] 

The initial sizing method was derived from Reference 12, 
which uses the average parameters of several air-to-surface 
missiles to determine the dimensions for the airframe under 
development. Once the initial sizing is completed, DATCOM 
(Ref. 20] is used to verify that adequate values for lift, 
pitching moment, and angle of attack are produced by the 
airframe. 

The munitions carrier design is driven by the restriction 
placed on the wing span parameter. From the desired wing span 


and equation 10, an overall length may be estimated. 


L =(L/b) Lb (10) 
= 95.2 inches 


The body diameter may be calculated from equation 11. 
= (11) 
D= L/(L/D)., 


= 3S] inches 
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Using the diameter, the avionics compartment and motor lengths 


can be estimated from equations 12 and 13. 


- [ ero 
<7 : reF/ 5 ) D il 
AVG 
= 32.8 inches 
gee es 
Louies ne a Dye wil 


= 27,14 inches 


Entering equation 14 with the payload weight requirement of 


500 pounds, the carriers gross weight may be computed. 


W, a Na/ dil cis) ee Oee 


ANG 
= 1150 lbs. 


The total lifting area is determined from the desired wing 


loading and gross weight of the weapon. Equation 15 is used 


for this determination. 
Sus We /(W/S) ave (15) 
= \0.3 f+! 


This area 1S a combination of the canard and wing surfaces as 


shown by equation 16. The individual areas are determined by 


dividing equation 16 through by the canard area. 


6. > : on (16) 


5.7 Ce ae 2 SiS lie 


47 


The overall length was determined from missile 
configurations that carry their payloads internally. In the 
case of the munitions carrier, the payload is external to the 
Main body. Therefore, the length can be adjusted to 
compensate for this difference as shown by equations 17 and 


18. 


exe 


Cy OP ed ee 


= Z\il.6 inches 


VG 


Ce wee ee (18) 


= 73.6 inches 


The perceived drag on the munitions carrier is assumed to be 
higher than that of a conventional design because of the 
carrier's external payload. The dimensions of the payload are 
nearly the same as that of the carrier. Hence, drag may 
nearly double that of a conventional missile. This increase 
in drag should be reflected by adjusting the length of the 
propulsion section. 

A Simple analysis shows that doubling the drag increases 
the motor length by 45% if the diameter is held constant. 
Therefore, this adjustment causes the overall length to 


increase to 84 inches. 
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The canard and wing dimensions are calculated next. A 
mono-wing layout was selected. This layout was chosen over 
a cruciform layout because the external payload may cause 
sufficient interference on the lower wing panels to effect 
their ability to produce lift efficiently. To meet the wing 
loading requirements, the wing planform must have a sizable 
area for a modest wing span. A clipped delta wing satisfies 
this necessity. Even though a delta wing is not as lift- 
efficient as a rectangular wing in the subsonic regime, it 
allows a wing folding mechanism to be avoided. 

To select a desired clipped delta wing, a taper ratio must 
be specified. The taper ratio,A, 1s the ratio of the tip 
chord to the root chord. For a given span, as this ration 
increases the wing lift curve slope decreases. A value of 
0.25 was selected for the taper ratio. 

pieesCou chord, CC’, and tip chord, C., may be determined 
from the wing span and area as shown below. Figure 9 defines 


the geometric parameters used. 


ay ft oe (19) 
ae CG Be 
as 4Cy 


Therefore, 


From Figure 9, 


Sw=0.5C,b+ &b a) 


where, 


a i a 


b/2 


Wing Geometry Definition 


Figure 9. 
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Rmecelecting agtip chogdlength,..a.root..chord,Jength may be 
determined from iterating equation 11 with the wing area equal 


to 9.5 ft*. Table 13 displays the results of this iteration. 
Hence, é - 1.4 ft 
¢t e 
5.4 ft 


Table 13. C. Determination Results 





The canard area was previously defined by equation 16 and 
an examination of Table 12 reveals the average canard aspect 
Bacion is 4.7. The planform selected for the canard 1s a 
delta wing with a taper ratio of zero. Consequently, the 
Canard span and root chord length may be determined from 


equations 21 and 22. 


b =v, AR, (21) 
24 inches 


ae i om Syl (22) 


= 10 inches 


aval 


With the parameters determined a line drawing may be 
drafted, as shown in Figure 10. As can be seen from Figure 
Oe the wing has ae Significant leading edge sweep 
(approximately 60°) with an aspect ration of 0.84. To enhance 
the lift curve slope of the wing, the aspect ratio should be 
increased. This increase may be achieved by reducing the 
lifting area by 20%, along with a reduction of the maximum 
gross weight allowable by 5%. These alterations increased the 
wing loading by 15%. The elevated wing loading is considered 
to be well within the structural limits of the wing. A new 


aspect ratio may be determined from the recalculated wing area 





as shown in equation 23. If, 
S$. Ogham = 7 Of e* (23) 
then, AR 12 
om 
= \.03 


Figure 11 displays a suitable baseline configuration from 
which the preliminary design may be developed. 

The wing surface can be optimized to produce maximum lift 
for. Minimum drag. Redmon [Ref. 12:p. 177} discusses a 
procedure titled "lift-drag (F) function" which produces such 
an optimization. The F function was developed for missile 
design in the case where a restricted span is given and the 
root chord must be elongated to accommodate a desired lifting 


surface. This elongation may produce high drag for the lift 
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Wing Area Reduction on Baseline Carrier 


Figure 11. 


54 


generated. The F function procedure requires the user to 
select a desired wing area. This area is then combined with 
a range of aspect ratios to produce a set of unique wing span 
and root chord combinations. The wing lift curve slope and 
zero lift drag are calculated for each aspect ratio. The F 
function, for each aspect ratio, is resolved from the 
summation of the normalized zero lift drag and the inverse of 
the normalized lift curve slope. The zero lift drag is 
normalized with the maximum zero lift drag value of the 
selected aspect ratios and ene lift curve slope is normalized 
with the minimum lift curve slope calculated. The F function 
Values are plotted against the normalized mean aerodynamic 
chord of each aspect ratio. The graph produced has a single 
minimum. The value of the F function for a given aspect ratio 
remains close to that of the minimum over a large range of 
aspect ratios. Therefore, selection of an aspect ratio within 
this range will produce a near minimum drag for a maximum lift 
Bemaition. [{Ref. 12:p. 178) 

To exemplify the above discussion, the following procedure 
is described. The F function is determined for the baseline 
derived wing surface area of 7.8 feet’ over an aspect ratio 
mange Of 0.2 to 3.0. The aspect ratio of 1.03 is highlighted 
for demonstration of the procedure. 


elven: 


Sr OE 
AR= 1.03 
A= 0.25 


aye) 


then, 


b,=WSVARY = 34 inches 


Ch= 25,,/Cb(itAl] = 52.9 inches 
€ = plage = 27 inches 


The calculations above are repeated for each aspect ratio. 
Upon completion of this step, the zero lift drag is then 
determined. DATCOM [Ref. 20:p. 4.1.5.1-2] provides an 
equation for estimating the wing zero-lift drag. Ri s 


equation is valid up to the transonic regime. 


O° & [I+ L(:}+100(£) | R (oe | 3 
where, 


— is the average section thickness to chord ratio 


\ec+- 


(CC 1S the skin tricti1on coc sete teme 
Ris the lifting surface correlation factor 


ie is the airfoil thickness location parameter 
The wetted area of the wing may be determined from 
equation 25 and S,, is the projected area of the wiwanag 


plan orn. 


8 et” 2(S..)— (Cy, D ) (25) 
= 12.6 fe? 
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Entering Figure 12 with the design Mach number, M, and 
wing sweep angle measured at (t/c),,, will produce R,,. The 
sweep angle may be determined from equations 26 through 32. 
The geometric parameters used in determining sweep angle are 


defined in Figure 13. 





Figure 12. Lifting Surface Correlation Factor [Ref. 20] 


GAC, (26) 
13.Z inches 


ae (27) 
Ge 


= 
fi 


i] 


Y= mile. ) (28) 
= 31,7 inches 


Si) 


ei aa) (29) 
= 7.9 inches 


Q= Y-Z (30) 


= 23. B3inches 
G=b/2 (31) 
= 17 inches 
UL. = ARCTAN (% } (32) 
= 544° 


From Figure 12 with an M of 0.5 a value of O33eu 
produced for R, « 

The C, may he determined by using a procedure shown in 
DATCOM [Ref. 20:p. 4.1.5.1-26]. This procedure requires the 
Reynolds number, R,, for the wing to be calculated andvemmeaes 
into Figure 14.0 Phe wk.) Vancom. 10’ based on MAC. This 
value must be compared with the Reynolds number cutoOi iia 
R,,, 18 determined by entering Figure 15 with M and cieuaias 
ratio, where 1 is the reference length and k is the surface 
roughness height. The value for k is determined from DATCOM 
[Ref. 20}. For this problem 1/k is equal to 9.2 & fo | 
Therefore, R,.. is equaliemer oa 10°. This value is indeed 
less than R, so it is used to enter Figure 14. This opemamaen 


yields a value of 0.0033 for Cy. 


a6 





Figure 13. Wing Swept Angle Geometry 


aye, 








Figure 14. Turbulent Mean Skin-Friction Coefficient on an 
Insulated Flat Plate [Ref. 20] 
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Figure 15. Cutoff Reynolds Number [Ref. 20] 
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With the maximum airfoil thickness ratio occurring at 
0.40, L will have a value of 1.2. This term and the thickness 
ratio are entered into Figure 16 to determine the pressure 
drag (PD) where, 


Po= [1+ uft}riooft} ] 


For this casesthe valuesrs.ieae 


1.2 for (t/c) located at x, 2 0.30c 


max 
OT (t/c) located at xs 0.30c 
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Figure 16. Subsonic Wing Minimum-Drag Factor [Ref. 20] 


Inserting the derived values into equation 24 yields, 


C, = 00056 
oO 
The wing lift curve slope may be estimated from DATCOM 
(Ref. 20:p. 4.1.3.2-1]. This estimation was derived foamea 
three dimensional subsonic wing. Equation 33 displays the 


formula: 
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2 2 
Caer oii G2: 
C = AR(2m)/[(2+| AR A = ) |+4)] (33) 
0% - 
where, 
a i | Oe M*) ” 
Ses, is determined from equations 26 through 32 
Substituting into equation 33 yields, 
(C, ) = 1.41/ radian 
a WwW 
The function may be determined from the above variables. 
The variables for each aspect ratio are listed in Table 14. 
From Table 14, the F function versus the normalized mean 
aerodynamic chord may be plotted. Figure 17 displays this 
Graph. 


Table 14. F Plot Results 
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3.00 |4.0x10° |0. 0034/0.0070/21.8°| 0.40 [3.30] 0.26 /1.00/0.078|1.078| 
2.50 4.0x100|0. .0034|0.0068/25.7°| 0.48 |2.93] 0.28 |0.97/0.089|1.059| 
2.00 |4.5x10°/0.0034/0.0067/31.0°| 0.60 {2.49} 0.32 |0.96/0.105|1.065 
1.50 |5.0x10°/0.0033|0.0064|38.9°} 0.81 ie 0.37 |0.91|0.132 ae 
1.03 |5.5x10° 0.0033|0.0056/49.4°| 1.17 |1.41] 0.44 |0.90/0.184|0.984 | 
0.50 |6.0x10°|0.0032|0.0042/67.4°| 2.40 |0.71] 6.63 |0.60]0.366/0.966| 
0.20 |1.0x107!0.0029]0.0032/80.5°| 6.00 {0.29} 1.00 |0.46/1.000|1.460 
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Examination of Figure 17 shows the selected wing planform 
to be extremely close to the minima. Thus, the planform has 


been optimized to produce minimum drag for the lift generated. 


12 03 O04 O05 06 Of @2 ee 


C/C max 





ree 


Figure 17. F Plot 

A lift analysis was conducted to determine if the wing 
paired with a fuselage could provide adequate lift at small 
angles of attack. 

From DATCOM [Ref. 20:p. 4.3.1.2-1], an approximate value 


for the wing-body lift curve slope may be determined. Equation 
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34 coupled with Figure 18 produce the desired results. 


Poe, a 
where, 


lene m@ierr Lawro coerricle;nt 
For this analysis, the canard iS assumed to produce zero 


Mec. From Figure 18: 


(ie OG 


thus, the wing-body 1ift curve slope is equal to 1.37/radian. 
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Figure 18. Lift Ratio [Ref. 20] 

To find the required angle of attack, the lift coefficient 
must first be determined from equation 35, where L is the lift 
produced. 

= L/q 3) (35) 


and given, 


llOO Ibs 
7.8 ft* 
G= 306 Ib/ ft" 


wD T 
UI t" 
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Substituting produces, 
C= 0.35 
L 
Assuming (C,,),, to be linear for the purposes joeayaam= 


analysis, then: 


C= (CL ee (36) 
Solving for @ yields, 


= 146° 
The angle required is not considered tolerable for the carrier 
design. 
There are three conceivable ways to reduce a; they are: 
Bec er ease ie, 
- Ineseacerey. 
* “a COmbination ef ene above 
If the weight of the weapon is reduced, then less lift is 
required. Examining equation 26 shows a reduction in lift 
will increase C,. The only component that lends itseisneame 
weight reduction that could significantly decrease Cane 
motor. This reduction, however, reduces the weapon's range. 
Because of the span limitation, an increaSe in surface area 
is not desired. The enlarged wing area would decreasSeme@mean 
would also degrade the lift curve slope value by decreasing 
the wing's aspect ratio. An increase in cruise Mach number 
is an option. The dynamic pressure would be increased which 
would produce a lower C,. An analysis was done over a range 


of Mach numbers to examine the effects Mach number has on a. 
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The results are shown in Table 15. Even with a cruise Mach 
number as high as 0.8, the angle of attack is still 
unacceptable. 


Table 15. Required Angle of Attack for Given Mach Number 





Therefore, it becomes apparent that the munitions carrier 
is unable to meet the operational requirements (OR) without 
increasing complexity (by adding a wing folding mechanism) or 
accepting a reduced range. 

Before proceeding with alterations that would allow the 
munitions carrier to fulfill the OR, the conventional missile 


configuration shall be examined. 


D. CONVENTIONAL MISSILE ANALYSIS 


Historical sizing was used to develop the conventional 
missile baseline. The baseline missile is determined by 
referring to Table 12 and recalling that an internal warhead 


weighing 150-pounds produces an acceptable kill probability. 
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The W, may be determined from equation 5: 
Wee (We/W \W 
a? (We wit de WH 
= A7Z2 \Ibe 


The length of each component may be calculated by 


selecting a body diameter, D. The diameter chosen is 11.8 


inches. 
Kay * (Lav/ Pye, P 
= 34 \nchés 
= 33 inches 


The warhead length 1s estimated to be 23 inches. 


Combining the results of equations 12 and 13 with the warhead 
length yields: 
2 scam en +L (a7) 
"ve oe M 
= 9Oinches 


The lifting area 1s determined from equation 15. 


SSN / ONE ) 
4.3446 


Equation 16 then allows the determination of the wing and 


\! 


surface areas, 
= + a 


=) = Sw 

a | + = 
then, 

5. = 0.34 Ft" 

S.= 3.964% 
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With the dimensions resolved, an analysis of the required 
C, can be conducted. Table 16 displays the derived lift 
coefficients for a selected range of Mach numbers. DATCOM 
[Ref. 20:p. 4.3.1.2-3] provides a method to estimate the value 
Seeche wing ybedy 1ift curve slope, (C,.),,- 


Table 16. Required C, for Given Mach Number 





(Cu) CKtK ot Kad (CL. lila) (38) 


alomeneeuarG otemose lift ferof the wing alone 
Kigy + Key, are the wing-body interference factors 
(ee See omc Iviom=s Ope Of the exposed wing 
Se opels) Che yratlolot exposed wing area to that of 
the total wing area 

The lift curve slope for the wing will be based on an area 

of 4 ft* and a cruise Mach number CieOn-S-. 
The planform selected is a clipped delta wing with a taper 
ratio of 0.25. From the above dimensions, the remaining wing 


parameters may be calculated. The results are listed in Table 


a7. 


OS, 


Table 17. Clipped Delta Wing Parameters 





These parameters are used to define the variables in equation 


Bree 
Solving for the wing-body interference factors: 
K K ( > \* (39) 
+ a 1+ 2 - 
WB) Bw) b 
= \,82 


K, Waly De det emiline cameo. 


Wee Giols ofr: (40) 


: (Cae Se 


where, 


CE is the lift-curve slope of the nose, 2/radian 


tad n 


and 


5 = wer 


= 0.79 ft* 
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PmceiiekN GeinbOwecuatLona3 94. yields: 


K= 0.33 


Therefore, 


(Cc, ) - Z.56/radian 


Using the above value and referring to Table 16, the angle 
of attack may be calculated. Entering equation 36 with the 
appropriate values produces a = 2.9°. This a is considered 
satisfactory to Warrant further development of the 


Semi iguration. 


E. CANARD SIZING 


The procedure used to size the canard was taken from 
Roskam [Ref. 21l:p. 259}. This procedure varies the area of 
the canard to allow determination of the location of the 
missile aerodynamic center (XAC). The position of the missile 
center of gravity (XCG) is calculated and compared with the 
SAC, and the comparison shown in graphical form. The 
difference between XAC and XCG produces the static margin. 
Entering the graphical comparison with a desired static margin 
yields a unique canard area. Selecting a static margin equal 
to 5% of the overall missile length will reduce the work load 
required from the control system because the missile becomes 
inherently stable [Ref. 19:p. 94]. 

Construction of the graph used for canard sizing is 


developed by initially determining the xXCG. Component weight 


Gk 


determination may be drawn from weight estimation techniques 
outlined by Redmon [{Ref. 12:p. 258]. Assuming the avionics 
system is to be less complex than those for the missiles 
listed in Table 12, the volume of this section may be reduced. 


This reduction reduces the length of the avionics section by 


10%. 
5 ie Guidance/Seeker Weight 
Wee” 0.00485 (W,) (Leg) 9 (41) 
where, 
L. 18 the dengtheorsenesceceia. 
De. 1s the idlanet creo Raine. Sermon 
eMekie 


Wass rain eyes 


Using a BASIC program developed by Rabang [{Ref. 22] 
for motor sizing shows the motor length may be reduced by 8 
inches. The reason for the reduction is that the thrust 
requirement is estimated to be lower than those for the 
missiles listed in Table 12 because these missiles are 
required to either fly longer ranges or cruise in the 
transonic regime. 

Appending the missile with the above adjustments 
yields a new length of 79 inches. The remainder of the 
component weights may be calculated based on the adjusted 


length. 


v2 


ae Weight of Body shell 


0.64 halt 
W-= 0.0604(L-) (D-) (42) 
where, 
L, is the body length 
Deetseehne body diameter 
W.= 89.5 Ibs 
oe Weight of Single Wing Panel 
|.02 0.S6 
Ww = 77483 (Sp ) (AR) (43) 
where, 


Srseene Ene e pare exposed surface area 
Res She <engite panel exposed aspect ratio 
Wy = lO lbes 

Since the wings must be detachable for storage, extra weight 
must be added to account for wing mounting brackets. These 
brackets would not be needed if the missile were of a uni-body 
construction, however. 

Table 18 displays the weight, moment arm as measured 


from the nose apex, and moment of each component. XCG may 


then be derived from equation 44. 


KCG= 7 WE (44) 
= 49 inches 


The XCG value is divided by the wing mean aerodynamic 
chord (MAC). This value is then plotted against the selected 


range of canard values. Because the canard weight is minimal 
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amo | |__| 


when compared to the gross weight, the XCG 1s assumed to be 





unaffected by increasing the canard size. 
To determine the position of the missile aerodynamic 


center, equation 45 is used. 


KEK, 0 (SEK Weill F cs) 


wR c 
where, 
F = 
Se eo = @ 
Ce ss | om, sa V7 Ue a y - Sel (46) 
The lift curve slope for the canard may be computed 
if a planform has been selected. For this case, a delta 


planform was selected. An aspect ratio (AR) of 4 with a taper 
ratio of zero is selected for the canard. Entering equation 
24 with the above parameters yields a canard lift curve slope 


value of 4.4/radian. 
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The aerodynamic center of the canard may be 
determined from Roskam [Ref. 23:p. 306]. Entering Figure 19 


i. 
a AR TANA,, = 4 


B/TANSA ip = 0-6 


a position for the aerodynamic center (XAC) may be extracted. 
The output multiplied by root chord gives the XAC measured aft 
from the canard apex. FOme this gease, XAC/C. = 10.6. The 
distance selected for the canard XAC when measured from the 
reference point (RP) is 41 inches. The RP is defined as the 
leading edge of the wing MAC. The distance of 41 inches was 
selected to move the canard as far forward as possible without 
impinging upon the nose cone. For use in equation 36, the 
distance is divided by the wing MAC. 

Next the wing-body lift curve slope is determined. 
Equation 29 is used for this determination. The nose lift 
meve Slope used to caiculate K, was estimated to be 
2.0/radian. This estimate is based on slender body theory 
which assumes an ogive nose. For this missile however, a 
blunted cone is selected, as it is less expensive to 
Manufacture than the ogive. The trade-off is that the cone 
is less efficient at producing lift [Ref. 13:p. 277]. DATCOM 
(Ref. 20:p. 4.2.1.1-1]) was used to refine the nose lift curve 
Slope value. To use DATCOM a fineness ration must be 


determined. The fineness ratio is defined as the length of 


1D 





Figure 19. Effect of Aspect Ratio, Sweep Angle and Taper Ratio 
on Wing Aerodynamic Center [Ref. 20] 


the nose to the base diameter. Assuming an infrared seeker 
will be used, Lindsey [Ref. 13:p. 278] recommends a fineness 
ration of one. See Figure 20 for nose cone details. 

For a cruise M of 0.8, DATCOM requires calculating 
the lift curve slope of the nose for both subsonic and 
supersonic flight regimes. These values are faired to produce 
the lift curve slope for the transonic Mach number selected. 
The value calculated using this procedure was 1.8/radian. 


Inserting this value into the equation for —K, andmmmen 
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Figure 20. Blunted Nose Cone with fn = 1 
calculating the wing-body lift curve slope yields a new slope 
fee. 53/radian. 

DATCOM [Ref. 20:p. 4.3.2.2-1] was used to determine 
the location of the wing-body aerodynamic center (XAC). 


Equation 47 is used for the determination. 
Xpae ( ‘ic y ; é 
Xac ( TAC: N aan "re wip) gee Loe wi8) KB) @ 


bre Sonar + Ce wins 7 CLE Biw) 





See Figure 21 for parameter definition. The procedures used 
for determining the values of the factors in equation 47 are 


shown below. 


(ve = Fae uy: (s) (48) 


oy 


tel. 


(Cy) eww K ow Cen de (Se/$ w) ~’ 
(Ciy ly = CH Gabe? oer (50) 
where, | 
(CCS a7 Sue 
Ki») and Ki, ane determinedmt rongh ign ne 22 amar 
= |,30 
Kw (B) 
K gi 0-52 
Therefore 
= |,55 
CLedvig 
(CL eran 0.63 


Summing the results from equations 48, 49, and 50 produces a 
Value of 2:53 pemTadmangtor ce 


Figure 23 is used to determine 


a 
Cre / we) 





Wate hls 
AR= Zz 
TANAL ine: 
»~ = 0.25 
yielding 
[Yas 0.55 
Ce 
e€ 
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exposed wing area 





Figure 21. XAC Geometry Definition [Ref. 20] 


To calculate 





Ae | 
Ch, | BW) 
Smemvalue Of the term, A.B, must be resolved. If this term 


is less than four, (Xa | Pemecdhemraced £Omn A.§ = 4 and 
Cre Biw) 
AB = 0, then interpolation is used to derive the correct 


value. For this case, 


AgBe 1.0 


qs) 


1.8 


r CEES ee = 
pt 4 HH 
atti ttt | Trt [ | | dT 

Ea | by) | 





— + it = t tt - 
y 
| | | 
ore pe - 
Hy eitea al 


body diameter 


me ae Se 


wing span 


a 
b 





Figure 22. Lift Raeres Kuice) and Key) 
Fixed Incidence--All Speeds. [Ref. 20] 


--Slender-Body Theory-- 
Calculatings ior A> aan 
X =e oe 
Cees) et (Zee TNAy ) 
re Bw), Te 
where G is the wing-lift carry over on the body parameter. 


This parameter is selected from Figure 24. Therefore, 
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ama if 





eo) 
ee AS 
Nong 4 
then substitution yields: 
ag) 2a SS) 


Dome ene mecdoce Olen pe — 0), Figure 25 15’ entered with 


W= a [ Acl!+Ae) TANA, ¢] 


then, 
W= 0.66 
migure 25 yields: Y 
[Moe = 0,33 
Cre Biw), 
The interpolation produces: 


=0,35 





| es) 


ine fanal term to be calculated is 


(el 


where, 
teal oo [se ( Lea.) (51) 
Cre N feq N Cre 

The first term is extracted from Figure 26. With a nose 


fineness ratio of one, 


(Kee) = -0,33 
oe a —_ 
The second term of equation 51 1s calculated by dividing the 


nose equivalent length by the exposed wing root chord. The 


nose equivalent length is defined as: 


fey, = 96a D 


ade 










r= 
Sar See SAM 
_| | | | | ae [| — 
Be? ) ieee 
Powe ae eae 







Figure 23. XAC/C.. Graph [Ref. 20] 


where, 
Fed = Lael 6 F gepy 
= 7 
Therefore, 
Keq= 84 
which produces, 
(Fe4) = 4.2 
ane 
Multiplying 
(7 | 4 
Leg on 
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yields, 
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Figure 24. Parameter Used in Accounting for Wing-Lift 
Carryover on the Body [Ref. 20] 


Substituting the above values into equation 47 bears, 


(tas\ wo,22 


Cre WB) 
Thus, the wing-body XAC lies 3.2 inches aft of RP. 


For use in equation 45, this length is divided by the MAC. 
The final value to be computed for use in equation 


36 is the upwash gradient, dé - Roskam [Ref. 20:p. 274] 


dm 
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Figure 25. Theoretical Aerodynamic-Center Locations for BA 
[Ref. 20] 


resolves this parameter graphically. The graph is displayed 
in Figure 27. Entering Figure 27 with the distance the canard 
lies forward of the quarter-chord point of the wing root chord 


in units of rootwehord . ytelde. 


de 


— =0.05 
dx 


The canard surface area will be varied from 0.5 to 
1.0 ft*. The combination of this range with equations 36 and 
37 produces Table 19. These values are plotted to generate 
Figure 28. The difference Between XAC for the missile and XCG 
when multiplied by MAC produce the static margin. If the 
static margin at launch is selected to equal 5% of the total 


missile length, the missile may become too stable as the 
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Figure 26. Aerodynamic-Center Locations of Various Noses 
(Slender=-Body Theory) [Ref. 20] 


flight progresses. This increased stability is caused by the 
reduced motor weight as the propellant is burnt off [Ref. 
inp. 302]. Therefore to avoid this problem the launch static 
margin selected is 2% of the overall length. This selection 
ensures the missile maneuverability towards the latter phase 
feet light. The gross canard area that equates to the 2% 


static margin requirement is 0.734 ft’. 
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DISTANCE FORWARD OF ROOT QUARTERCHORD POINT IN ROOT CHORDS 





Figure 27. Wing Upwash Gradient [Ref. 20] 


Table 19. XAC for Selected Canard Area 
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Figure 28. X-Plot Graph 


F. DETERMINATION OF THE WING-BODY-CANARD LIFT CURVE SLOPE 


The lift curve slope for the missile was examined to 
determine if the cruise angle of attack would fall within 
acceptable limits. DATCOM [Ref. 20:p.4.5.1.1-2] furnishes 


equation 52 for such a determination. 
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: er (S\ (5 S28 
3) 21), [Ket Kenko ie! 


€ ao We 





(Ly = (Ct anes K staph 


This equation is valid for: 
+ M up toe and” ine ludingee-c 
+ Wing span to canard span ratio greater than 1.5 


+ Body diameter to wing semi-span ratio less than or equal 
iO e028 


The third requirement is not met by the canard, as this ratio 
is equal to 1.1. Therefore, the calculation will be completed 
with the understanding that the canard area needs to be 
modified to meet this requirement on the next iteration. 
Perhaps, to avoid this situation during future design efforts, 
requirement three could be applied to RoskKam's canard sizing 
method. If the requirement is not met, then the canard area 
derived is the exposed area versus the gross area. 
Examination of the first term in equation 52 reveals that 
the interference factors must be based on the canard area 
instead of the wing area as previously calculated. Replacing 
the appropriate terms in equation 31 with those that are 
relevant to the Canard) 9y lelds eee eooe The remaining 
interference factors are determined from Figure 22 entering 


with a diameter to canard span ration of 0.58. 


Mey 


K Biw)- 0.94% 
Therefore, 


CK +K |= 4. 2G 


wey” een 
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All variables in the second term of equation 52 were 
previously determined except for the dynamic pressure ratio, 
ol Ae! For early stages of preliminary design, the dynamic 
pressure ration may be set to unity (Ref. 23:p. 274]. 

Substituting the derived values into equation 52, 


produces: 


C, = le-27radian 


O¢ 


2 Orac/ 


Piemmectined —Cg@ibased on the canard area is found from 


equation 35. 


es 


CL eb /4 Se 
=0.69 
From equation 36, the unknown angle of attack may be derived. 


Rearranging, 


Ce Celera 


CR 


SE 


The increase in drag due to this angle of attack should 
be minimal; therefore, the lift curve slope value is 
considered acceptable. 

The question of stability must be addressed to ensure an 
adequate control system may be realized. The next section 


concentrates on this area. 
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G. DETERMINATION OF THE PITCHING MOMENT CURVE SLOPE 


This section will determine the longitudinal static 
stability of the missile. Anderson [Ref. 24:p. 354] states 
the requirements needed for longitudinal stability (LS) are 
ena cr 


- ¢, MUSt De wegcenive 


- Cy, must be positive 


where C,, is the pitching moment curve slope and C)Miieueege 


Value Of CO Sate Zor cm Jace Cc, 1s defined as the pitching 


M 


moment coefficient. 
DATCOM [Ref. 20:p. 4.5.2.1-1] provides equation Sammie. 


cleheeneneloibiale; ()- 


aes fap oly , acy a { 5 m : a 4 \, a 
O, hae eae =) A, Tae K ewe Clee t \s, \(= (53) 
Gs / 


x 


where, 
A= Oo (SH (Se) ~ Se - 0.28 ee 
Le eg! a Ce 
and 
A > es XG 
fs a es (55) 
- ox ore 
Bugure 22> derines G7 sand ale 


The term (Ky / Cy. ) is defined by equation 56. 


(56) 








Nee = ( Xe | Cre, i D TAN 
\ 4 « = LE 
Cie Ten Ciae Ly 
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where (Xec / Op ) is determined from equation 38 with the 
es 
Canard area as the reference area. 


fmus from Figure 29, 


oe = 0.44 


arc.) = 12:57 

(9./ a= 6.25 

Oey G.) =- 5.14 
Substituting into equations 54 and 55 yields: 

A oa wa 

dee = oO. 74 


é 


Entering equation 53 with the above determined values with 


those values previously calculated bears: 


6 -5.9fradian 


Mo 
= =o) © 


The canard area is the reference for oan ig Cy. is 


desired with regards to the wing area, then the coefficient 


is multiplied by the canard area to wing area ratio. 


Ca, ($2) = -1.1 /radian 





The first LS requirement is met; C, may be investigated. 


C,, can be considered the y-intercept on a graph of C, verses 


Mo 


a. Therefore, 


C 


= + C fh 575) 
P C4 Ae 


out 





Figure 29. Vehicle Geometric Parameters [Ref. 20] 


a2 


ema to determine C,, C, is set to zero. Thus equation 57 


reduces to: 


Cy 


)) 

GS 
S 

Pas 


.) 
where, 
<a. 2.4: 
then, 
Cy, = 0.05 
The results of equation 57 meet the second longitudinal 
stability requirement. Therefore, the missile is 
longitudinally stable. 
itiemscanard incidence angle, af, needed to develop the 
Beetoed Moment COerricrent, C., necessary to hold the missile 
at the cruise angle of attack will be determined. The wing 
is assumed to have a symmetric cross-section with the 
Mmeeaence angle, i,, set to zero. Roskam [Ref. 23:p. 320) 
eeipelies equation 58 for determining ¢C,.. 
ae Pls Crem (58) 
where, 
eels ene Zero-litte moment of the wing-body 
combination 
(ey ace eneeze~o-litt moment of the canard 
When the wing-body combination is at the zero-lift condition, 
it will be assumed that the moment created is negligible. 
Thus equation 58 reduces to: 


Cue (C A 2) 


Cc 


The zero-lift moment of the canard is defined in equation 59. 


\ \ Rees 
(Cade = (oe eee (59) 
where, 
(C ) = C ami t& ) 60 
Bole Le Sw C ° (60) 
The canard down-wash angle, €, , is assumed to be zero. 


Substituting in the variable calculated earlier, bears: 


C.). = Ca (62) 
Recalling that, 


Rie, = ace 


then equation 58 becomes: 


Cs =(1b0)(0.81)t, 


0 


Solving for i, yields, 


C 


H. HORIZONTAL TAIL-SIZING AS DETERMINED BY STATIC DIRECTIONAL 
STABILITY REQUIREMENTS 


Roskam {Ref. 21l:p. 265] develops a procedure from which 
the horizontal tail may be sized as related to the desired 
yawing moment curve slope. The procedure uses a range of tail 


areas to generate a C,, value. Equation 62 is the governing 


Cn Coglet CL, canal AY) (62) 


formula used. 
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The geometric quantities are defined in Figure 30. DATCOM 
(Ref. 20:p. 5.2.3.1-1] may be employed to determine the first 


term in equation 62, where, 


26 \ en. / 
ae < S 
The wing-body interface, K,, is determined from Figure 31. 


Figure 32 defines the parameters needed to enter Figure 31. 


The values of these parameters are: 


Re = 6.6 ft 
Xm = 4.1 ft 

ae =) OFC Ie 
WW = h, = ae | 





Figure 30. Xv Geometry [Ref. 20] 
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Figure 32. K, Geometry [Ref. 20] 
Entering Figure 31 with the above values, yields, 
= 0.002 
Ky 
mie heynolcS number CTactor, K,,, 18 extracted from Figure 
33. The fuselage Reynolds number 
6 
R = 3.5x I0 
n 
yields, 
Keo we 
Substituting into equation 63 produces, 


(Cope SOLOS. 


To calculate the lift curve slope of the horizontal tail, 


‘ent maeene and 2/\.. wiliweee estimated since once an area 


av! C/2 
is selected, that area may be molded to fit the estimated 


Pemameters..» If, 


a7 


=e i a 


a iv, | aaa eee —— = 
+ 1° 
* 





wo. 
oO r=) Oo rid 
eg m ON at 
mn 2 
YY =-_—* 
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Figure 33. 


Effect of Fuselage Reynolds Number on Wing-Body Cu 


20] 


[Ref. 
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it 


1. G2 


Or 


AR 


C/2 


oy. = 2.43/radian 


!) 


then, 


The moment arm of the vertical tail is estimated at 20.5 
inches. Equation 62 becomes, 

Cup = -0.009 +0.0064 5, (64) 

Roskam [Ref. 21:p. 266] recommends C,, = 0.001 for inherent 

directional stability. Substituting into equation 64 produces 


a tail area of 1.5 ft’. 


; DRAG DETERMINATION 


Throughout the design process, an attempt was made to 


increase the value of the critical Mach number (M The 


pe 
direction of this effort was to allow the missile to operate 
at the edge of the transonic regime where wave drag is 
minimal. Thus, the motor weighs less when compared to the 
motor needed to operate in a wave drag environment. All wing 
planforms have been swept to increase the critical Mach 
number; nevertheless, the nose may become sonic prior to the 
planforms. Thus, Raymer [Ref. 25] provides an estimate 


technique for the onset of wave drag. 


nt, 


2(f,) = 2 


29 


where f, is the fineness ration of the nose, then Raymer (Ref. 


25] gives, 
Me = 0.79 


Therefore a cruise M of 0.8 was selected. 
The drag coefficient for the missile in a cruise profile 
can be determined from DATCOM [Ref. 20:p. 4.5.3.]=3) an 


drag coefficient for the wing-body-canard configuration is 


given by equation 65. 


= 65 
Oa (Co) et (Co), + (Co) ot Ne ( ) 
where, 
€,, drag coefiicient  atiZeneSiieaE 
C,,; induced drag of the wing-body combination 


ACp= Cy, (Se/% Sy) 

Co= (Cp). + Gor lin 
The equation is stated to be limited to a Mach number of 0.7, 
but will be assumed to be valid for this design. 

Zero-lift drag coefficients may be calculated with the 
method outlined in DATCOM (Ref. 20:p. 4.1.5.1-2]. The initial 
step of the method calculates the skin friction coefficients 
of each component at M = 0.8. Figure 14 is entered with the 
Reynolds number, R,, for a given component. 

R, =7-V (reference length) /p (66) 
where, | 
V is the. velocuc, 


(? is the density 


100 


u is the absolute viscosity 


mnen, 


Ry, = 1.32 x10" 
RX, =  lnniaeee ip? 
a 5°22 2 New 
ey 


z 
Lge = Jade“ ic 


Figure 14 yields, 


6 = 0.0030 
fe 
= 0,909036 
Cr Zs 
a = 0.0025 
fv 
G =Meeco 25 
¥@ 


The component pressure drag, PD, is determined from Figure 
17 when entered with, 
Weenies 7 Co — O0.05,05 = 1.2 
Canam, Cc = 0.03, L = 1.2 


wa Mies Zeu— Ol, OS sels ele? 


paaoauces, 
ae CSC«OS 
PD. See 
PD, = 1.02 
The PD for the body is determined from equation 67. 
(Px) 
ie = { | + — , + 0.0025 4s | (67) 
= (Ne/d) o 
which yields, 
PD, = 1.22 


Je(8) 


Determining the lift-surface correction factor, RR) Seem 
Figure 12 bears the results in Table 20. 


Table 20. R Values for Components 


L.S. 





Next the wing-body interference correlation factor is 


drawn from Figure 34 when entered with R 


Rig 2 0-99 


The calculation of the wing-body zero-lift drag exclusive of 


ip thus genéeracumes 


the base drag, X, 1s considered in equation 68. 


[PD R= + SPDs (5) }Rup 8 


where, 


Cc, is the base drag coefficient 


Db 


S. 1s the body base area, 0.44 ee 


B 
S, 1s the surfacejarea jerebody,, ccm ft° 
Substituting the appropriate values into equation 68 gives, 


X= 0.0204 (69) 


The final term needed to complete equation 65 is Co, Which 


Cp 0.02.9 (78 S109 dp (70) 


1Sso.defined as: 


EOZ 





Figure 34. Wing-Body Interference Correlation Factor--Subsonic 
Speed [Ref. 20] 


where, 
Dis the base diameter, 9 inches 
(Co) is (Ce), (PD), (S5/ SQ) 
Jha 
ut) 
(Cp), \44. 


then equation 70 yields, 


Cp, = 0.032 


103 


The values of X and C,, are summed to produce (Co) oe as 


shown in equation 71. 
S 
(Gp) K+ Coy | Ys} (71) 
= 0.024 


The remaining zero-life drag coefficients are determined 


from equation s72, 


C= Cy (PD) R., (Se/ Sw) (72) 


substituting the proper values gives, 


(Code = 06x 10° 


(Cy,Jy = 20% 10” 


The induced drag coefficient for the wing-body combination 


may be calculated from equation 73. 
(Core7 Bot BY/AR+ BAR + B(TANA,.) + By (t/Z) + B.(f./p) (73 
+ eae BLA ae Be 7 oF Be A? + oe (Tp)+ B,(Lee/ C ) 
7 Doo os SCR Ys Coaleree ta Cu + Sys Sw 


) 


where, 
yo 0 oo are regression coefficients as function of 
Mach number and angle of attack obtained 
from Table Zi tOme ac. 
ope fineness ratio for areas shown in Figure 35 
Tp transition strip factor (=0 for this design) 
(Lens c ) ratio of the leading edge radius to the MAC, 
(estimated at 0.004) 
§ wing twist, (=0 for this design) 
Oana camber ratio, (=0 for this design) 
CL, conical camber design lift coefficient (=0) 
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Table 21. Regression Coefficients (B, - B 


On nw AW A) — 


OO DIM pm e&wDW — 


eS ee eee 
oOow~ mom BN — 


0 00221 
0 0033B 
001022 
0 02014 
0 03137 
0 04483 
0 05954 
0 OBS86 
0.12288 
0 18902 
0.19087 
0.19993 
0 21762 
0 22B25 
0.26345 
0 2672B 
0.26672 
0 38780 


B. 

-0 00011 
-0 000s9 
-0 00081 
-0 00136 
-0 00160 
—O0 00768 
—0 00275 
—0.00385 
—0 00533 
—0 00€51 
—0 00645 
—0 00359 
—0 00071 

0 0008 2 

0 00251 

0 00€15 

0 00963 
—0 01148 


Bie 
0 03464 
0 0£ 100 
0 07937 
0 10986 
0 12995 
017974 
O csa87 
C 269289 
0 34386 
0 4021B 
043511 
0 44284 
0 39269 
0 24896 
0 23602 
C 21889 
O.17675 
0 35079 


-0 00144 
—0.00302 
—0 00891 
-0 01193 
-0.01B95 
-0 02209 
—0 02427 
—0 04047 
-0 06032 
—0.094B4 
—0.12050 
-0.10093 
-0.05BBB 
—0.06174 
-0.08437 
-0 04519 

0 03024 
—0 06958 


B, 


0 00404 
001149 
001741 
0 02455 
0 03701 
0 04932 
0 06720 
0 09758 
0 13366 
0.14191 
0 12269 
0 12092 
0 12946 
01411291 
012110 
0 17305 
0.20127 
0 23079 


By 
0 03796 
0 05301 
0 08300 
0 11608 
0 64853 
0 01671 
—O 06230 
—0 20397 
—© 16500 
-—0 06673 
001019 
001175 
—0 55132 
003958 
-0 13912 
—0 40147 
—0 44420 
0 06490 


—0 00009 
0 00036 
0 00044 
0.00101 
0.00175 
0.00408 
0.00622 
0.00675 
0.00656 
0.00500 
0.00370 
0.00286 
0.00378 
0 00297 

-0 00007 
0.00045 
0.00162 

—0.00033 


B, 
—0.00852 
—0.02119 
~-0.03594 
-0 05123 
—0 08555 
—0 10395 
—O 13865 
0.1944) 
—O 27088 
—0.25370 
—O 19227 
—0.17696 
—0.23012 
—0.195B65 
0.33286 
-0 52066 
-0 66993 
—0 B 1556 


Be 
—0 00568 
-0 01138 
—0.01716 
—0.02432 
-0 03608 
—0 05082 
-0 06341 
—0 08062 
—O0 08937 
—0 09478 
—0.09377 
-—0 08753 
—0 08675 
—O0 07231 


_ =0 07000 


—0 08010 
—0 07856 
-0 07221 


-0.00073 
—0.00054 
—0 00118 
—0 00320 
—0 00608 
-0 00776 
—0 01036 
-—0 01503 
-001914 
—0.017B6 
—0.01996 
—0.009 38 
-0 01095 
-0.00059 

0.01070 

0.01068 

0 00616 
—0.04003 


By 
0 00366 
001016 
0.01916 
0 C2690 
0 04712 
0.05655 
0 07161 
0.09494 
0 13617 
011623 
0 0765B 
0.07638 
0 11647 
0 10526 
0 23396 
0.36192 
0 46675 
0 53140 


yg % 10 


OLS 


15) 


—0 00757 
—O 04293 
-0 09171 
—0 16796 
—0 76045 
—0 33466 
—O 46785 
~0 53205 
-0 61233 
—0 67543 
—0 70864 
—0.80847 
-0 87035 
—0.75253 
-0.72073 
-0B1612 
—0.73124 
—0.00836 


Bi 
=—0 00020 
—0.00084 
—0.00120 
—0 00164 
-0 00185 
—0 00274 
—0 00261 
-0.00314 
—O0 00244 
—0 00763 
—O0 00338 

0 00149 
0.00542 
0.00630 
0 01239 
0 01483 
0 C2226 
001sS3s 


Substituting into equation 73 produces, 


( CO ve = 0.00 ae) 
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BS 
0.00001 
0 00016 
0.00029 
9 00031 
0 00073 
0.00035 

—0 00004 
0.00026 
-0 00006 
0 00000 
0.00052 
-0 000B0 
~0 00312 
-—0 002B6 
-0.00351 
-0 00314 
—0.00435 
0.00680 


a5 
0 04956 
0.27954 
0 48014 
0 BB190 
1 46725 
22427 2 
252398 
Z2/6c2B 
3.34227 
411935 
4 66755 
3 29950 
441136 
2 .9ad00 
1.94911 
2 44009 
1 45434 
—4 53597 





The final term, C,, is composed of the Canard zeroa—imere 


drag and induced drag coefficients. The zero-lift drag 


————— 


coefficient was determined earlier and (C,.). 1s defined in 


equation 74, 


Z 
= (CG 
(ovo oie ( Le) /( fre AR) (74) 
where e is the Oswald efficiency factor, estimated from DATCOM 


[Ref. 20:p. 4.5.3.2-2} to be 0.5. The lift coefficienmiee 


defined as, 


OL (75) 


CL. = CK, 


Cc LKe 
where, 


Ke = OC gt ve 


Thus equation 74 produces, 


Introducing these values into equation 74 produces, 


( Uyeiee 0.0134 





Figure 35. l 1 


Definitions [Ref. 20] 


Ld A 


ING Ts 


Therefore, 
= 0.00067 0.0134 


= 9.0\4 


hence, 


AC» a Cy. ¢ on Sw) 
= (0.0|\4)(0,1864) 
= 0,012 GC 


Equation 65 produces a cruise drag coefficient of: 


C, = 0.034 
The cruise drag may be resolved from: 
D= 40,5. 
= (949)( 0.034) (4) 
= \20 \bds 


J. DETERMINATION OF THE PITCH AND YAW MOMENTS OF INERTIA 


The moments of inertia becomes important during the design 
of a control system for the missile. DATCOM [Ref. 20:p. 8.2- 
1] provides a procedure for estimating the moments of inertia. 
The assumptions made when using this procedure include: 
* two planes of symmetry 
- the CG lies between ( £/3) and (2573) 
* each component is solid (non-liquid) 
Figure 36 defines geometric parameters. The pitch and yaw 


moments of inertia (B,C) for each component may be determined 


207 


from equdadtvon=7e: 


So: _ 2 
ae 4 4\,Mm mi 3 23ypb , 2mIN-bT | 
BzC = (R125 Se Pee UC ei. eae 42 ) 


2 ~.,MN* © 2OTN FeN* ee 
+(f2— $7 (PES EETN y 4 FEN yg 


(76) 





26 | 





l 


Figure 36. 1,, 1, Definition [Ref. 20] 





where, 


F= (1, - 1)) 


mF 
II 


component weight 


>< 
It 


locationgoef compenenemeG wat .c-. 
the nose apex 
Ww, = (W/F) (4-[6(x - 1))/F]) 


w, = (W/F) ([6(x - 1,)/F] - 2) 


mM = (Woo 
T=) (de 
FW aD (Wa) eae) 
N-= a, — 1) (a) — 9d) ye) 
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b 


Weeeel, (Cw - w,) / F) 
ec (ce dj Fr) 
The lifting surfaces were not included in the estimation. 
These surfaces make up less than 5% of the gross weight; hence 
their effect should not significantly alter the outcome. The 


Calculations for individual body components are shown in Table 


eZ. 
The results in Table 22 are inserted into equation 77. 
=C =P Bt LWe*- (LW)Ke (77) 
where, a 
X= LXW/(LW) 
then, 


B= C= (26369)4( 1170814) — (442)(44) 7 
= 357|b£t° 


Table 22. Moment of Inertia Determination 





Veale! ale - | | DS aaa 


| section =6p{ J | K | fF ! W x | wo | ow 


| 2 
Cin) [Cin) [Cindy | Clb) | Cin) S b |ctb ft ) 
Pr ovecmidel)| sce S0e tere | 15 | 2.4 | 2.4 | 41.9 | 
| WARHEAD SOmeosme cS silears (41.5) |9725)| 7-5 |. 64.2 


PERORULGION oS ui7omi co (197 | 66 | 7.6 | 7.6 | 91.0 | 
a a Ee ES es 


K. FINAL BASELINE CONFIGURATION 


The missile's dimensions are listed in Table 23. Figure 


37 displays a three-view diagram of the overall configuration. 


Oe 


Table 23. Summary of Major Component Dimensions 


| 

| CANARD TAIL | 
2 2 2 

| AREA (EXPOSED) 2.1 ft 0.13 ft Ina ft | 


SPAN 34" 20.6" Stee. 


45 


| INCIDENCE ANGLE 0 het 
ROOT CHORD 13.1" (ExPOseD) 

| FUSELAGE | 
DIAMETER fee! LENGTH 79" 


if 1, 


TOP 


a = 





— a 
| 
SIDE 


aie le) Shou 





Figure 37. 3-View of Baseline Missile 
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VI. PROPULSION 


A. INTRODUCTION 


This chapter determines the dimensions of the solid rocket 
motor to be incorporated into the missile design. The main 
design tool was a BASIC program titled LPROP developed by 
Rabang [Ref. 22]. The program is based on a propulsion sizing 
method discussed by Redmon [Ref. 12]. This program@eag 
optimize the chamber pressure of the motor to produce maximum 
thrust from a given set of size parameters. The method 
assumes a constant acceleration boost and a constant cruise 
ale levce, 

The design requirements or goals for the motor are: 


- to give the missile a low speed launch (LSL) capability 
(>200 kts) 


- to provide sufficient thrust that will enable the missile 
to obtain a range of 15 nm. 


The first goal may be accomplished with the use of an 
adequate boost phase. A booster will accelerate the missile 
to the desired cruise velocity over a small time interval. 
In the case of this missile design, the lift provided by the 
structure at low speeds will retard the missile freefall so 


that the booster will not have to provide high g 


le 


accelerations. This capability also allows the delivery 
platform to launch from low altitudes (<1000 ft.). 

The second goal is obtained by proper sizing of the 
sustainer section. Once the missile is boosted to the desired 
cruise velocity, the sustainer provides adequate thrust to 
meet the cruise drag encountered. 

A baseline motor will be developed prior to the use of 
LPROP. This baseline gives an acceptable reference that 


should approximate the results given by the program. 


B. BOOSTER SIZING 


The maximum axial acceleration selected is 5g. This limit 
ensures a lighter airframe than that of a missile subjected 
to 30g's, provided all other variables remained fixed. The 
procedure outlined is that developed by Redmon [{Ref. 12:p. 
e7 } 

The boost time may be estimated from: 

tL =(v2-yi/a a 
where, 
Velo veMescruise velocity, fps 
Vv, is the launch velocity, fps 
a is the acceleration, ft/s? 


Equation 78 yields: 


a ei 
(394- 338)/lb61 


=e 


\) 


TAs 


If the specific impulse, Isp, is assumed to be 210 s and 
the cruise drag is set to zero, then an initial propellant 
weight may be estimated by equation 79. 

W,/ We = EXP (AV/g1,,) (79) 
where, 
W/W, 1s the launch weight to empty weight 7a 
AV = (v2-V,) 
Substituting into equation 79 bears, 
W/W, = ExP( 894 - 336 )A(22.2X 2107) 
= le 


Given a launch weight of 480 lbs, then: 


We = 480°/1.| 
= 437 \|be 
Thus, 
WE = We-We 
oom fol 


The total impulse, I for the booster can be estimated 


t? 


fLOlne 
=, = Dep Wp (80) 
— (Zieiie= ee 
= 7980 lobe 
The thrust provided by the booster is then: 
= (ZfAS2 Noe 


m4 


Peeeter volume 1s determined from: 


Vol, = Wi 7 (pe 20Ny) (82) 
where 
NL is the volumetric packing factor, 0.85 
Pt is the propellant density, 0.065 lb / in? 
HeENCce, 


Vol_= ©87.7 in> 


Provided a body diameter, D, has been selected, the booster 


length may be calculated: 
ait (4VolL) 7 (* D*) (83) 


I 


1) 


@.| \nches 


C. SUSTAINER SIZING 


The sustainer is sized by using the procedure developed 
by Redmon [Ref. 12:p. 187]. This section is required to 
propel the missile to its maximum range once the boost phase 
iomconcluded. The size of the sustainer depends upon the 
range and speed of the cruise phase. These parameters 
Sem@erate a burn time, t.. 

tee CRS Wale (84) 
where, 
R is maximum range 
Peto ees taneommbamercdsduring the boost phase 


We SUES NS (eet ae hetero 


aS 


From Newtonian physics, for a constantly accelerating 
body, the distance traveled is: 
R, =WVa)at® (85) 
= (1/2) Wei 3.5)" 
= 9beF+t 


The maximum range required is 9.0 x 10° ft and the cruise 

velocity is 894 ft/s. Substituting into equation 84 yields: 
te= |006 

For the cruise phase, the thrust is equal to the drag, or 
T.= D, 

During the design of the airframe a drag of 130 lbs was 

calculated. This drag figure is increased by 25% to account 

for the possibility of under estimating pressure drag. 


The sustainer impulse can be calculated from equation 86. 
L.- Dt. (86) 

(\bS)(100) 

PIELe” = 


The propellant weight of the sustainer is determined from: 
We = s/ Lcp (87) 

(16900)/ (2° 

TY \bs 


Following the same procedure as used for the booster, the 


i\ 


)} 


tt 


volume and length are determined: 


Vol, = 1208 1m? 


L. = (oainene= 
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D. NOZZLE SIZING 


Redmon {Ref. 12:p. 207] provides a method for sizing the 


motor's nozzle. The nozzle throat area is determined from 


equation 88. 
are 
A= Yc, P) (88) 
where, 
Gres eneren msc coefficient 


Be isethe Chamber pressure 


The thrust coefficient is determined from equation 89. 








a rik 
emcees | 














eens 5 (89) 


C 


C.. = ye, 


_ oy 


| — 
where, 

ecm. hozete ert ecrency 

% is the nozzle half angle correction factor 

¥ is the specific heat ratio of the propellant 


P, is the ambient pressure 


0 
Equation 90 provides a solution to the nozzle half angle 
correction factor. See Figure 38 for a geometric definition. 
A= (1+ CosK)/2 (90) 
Redmon [Ref. 1l12:p. 192] recommends: 
@ = 15° because larger values give significant non-axial 
flow components and smaller values tend to create 


elongated nozzle lengths. 


Therefore, substituting into equation 90: 
A= (1 4 COS 156°)/ 2 


= 0.98 


aly 


nozzle throat 


nozzie exit 





Figure 38. a@ Geometric Definition [Ref. 12] 


The nozzle efficiency is shown from historical data to be 
about 0.96 [Ref. 12:p. 191]. The ambient pressure is taken 
at 15 psi and chamber pressure desired is 200 psi. The 
specific heat ratio for the propellant is selected from Table 
24. The value chosen for this caleulations1s 127 ee 

Substituting into equation 89 produces, 

Ces 1.94 


Thus equation 88 becomes, 
A, = 0.4 1m? 


The nozzle exit area may be determined from the ratio: 


A -L§ % Neo TET D7 2 (G1) 
7h, Bes | i+ ( 2 M i ‘4 
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Table 24. Propellant Characteristic Values [Ref. 12:p. 













a 


Therefore, 
Ae= 3.4 in 
The nozzle length is found from: 
ae = De - De 
2aAN See 


then, 
Ln= 2.6 inches 


198] 


(92) 


To optimize the chamber pressure, propellant weight and 


caSing thickness, the program was executed with the results 


shown in Figure 39. A cross-section of the final motor 


configuration is shown in Figure 40. The final motor design 


meets the initial design goals and can be packaged within the 


dimensions of the missile body. 
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Final Motor Parameters 


Figure 39. 
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Cross-Section of Motor 


Figure 40. 
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Vil. SEEKER DEVELOPMENT 


A. INTRODUCTION 


For an autonomous weapon to intercept a desired target, 
the weapon must have an effective homing guidance system. The 
term homing guidance is defined as "a missile guidance system 
that can determine the position of the target and can 
formulate its own commands to guide itself to the target" 
[Ref. 26: passnezip Formulation of the guidance commands is 
covered in the autopilot design chapter. This chapter will 
address how the missile can determine the target's position. 

Two types of homing were considered: they are active and 
passive. The active homing system would use radar mounted 
onboard the missile to locate and target the intended target. 
The passive system would use an infrared (IR) detector for 
target tracking. The pasSive system was preferred for this 
design, as it has two main advantages over the active system. 
The first advantage is that the missile would not transmit a 
Signal that may alert the intended target of an impending 
attack, thus reducing the delivery platform's susceptibility. 
The second advantage is that the infrared system is less 
complex, and would cost less than the active system. (Ref. 


262 >. LOSg 
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A definition of terms will proceed an explanation of the 
design process. These terms and their definitions are taken 
from Masters [Ref. 27:p. 2.2] 

Radiant Flux : P, measured in Watts, is the radiant energy 
per unit time emitted by the source. 

Radiant Emittance : W, units (Watts/cm’) is the radiant flux 
per unit area emitted by a source. 

Radiant Intensity : J, units (Watts/steradian) is the 
radiant flux per unit solid angle emitted 
Paomra Source. 

Radiance : N, units (Watts/cm*--steradian) is the radiant 
emittance per steradian emitted from a 
source. 

Lmeeaeiance > H, units (Watts/cm’) is the radiant flux per 
unit area incident upon a surface. 

Transmittance of Atmosphere A, 1s the ratio of the 
radiant flux transmitted through the 
atmosphere to the radiant flux incident to 


the detecting material. See Figure 41. 


B. APERTURE DETERMINATION 


For the detecting material to become exposed to the 
radiant thermal energy, the energy must first pass through the 


seeker's aperture. Therefore, the lower the radiant emittance 


2.3 
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Figure 41. Transmittance of the Atmosphere [Ref. 27:p. 10] 


from the target, the larger the aperture must be for a fixed 
detection threshold. This relationship is analogous to the 
human eye; as light intensity 1s decreased, the iris opens the 
aperture or pupil to permit more light to enter. However, in 
the case of a seeker, the mechanical actions of an iris would 
be costly; hence, for this design, the aperture is assumed to 
be a fixed area. 

The aperture area depends on the targets thermal 
Signature, the sensitivity of the detecting material, and the 


range from which the seeker views the target. 
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1. Target Signature Estimation 

The thermal signature is developed from a fictitious 
target with a hull of 500 m’ and a body temperature of 20°C 
(Ref. 28]. 

The target can radiate energy across several bands 
of the electromagnetic spectrum. The band at which the seeker 
is to search for the target must be resolved. To assist in 
the selection, Birk [Ref. 28:p. 2] states: 


The atmosphere absorbs IR radiation except in certain 


regions of the electromagnetic spectrum. These regions 
of the spectrum through which IR radiation can pass are 
known as atmospheric windows. Two important windows are 


located at approximately 3-5um and 8-12 um ranges. 

Gates (Ref. 29:p. 136] compares the trade-offs of 
each window. 

The 8-12um window renders IR decoys ineffective and 
makes it difficult for the target to camouflage itself. 
However, this window is susceptible to environmental 
hactors. 

Birk [Ref. 28:p. 3] adds that: 

Objects in the 8-12um tend to show the entire body and 
therefore the image is extended. .. .IR missiles tend to 
use the 3-5yum range so that they can home on hot spots 
which appear as point sources. 

Despite the advantage of having the target presented 
as a point source, IR decoys are very effective against such 
seekers because the decoy presents itself as a point source. 
Peer. 29:p. 136]. 


The 8 to 12um window is selected for this seeker. 


The detecting material will have a geometry that will not rely 


25 


on point targets for homing information. Therefore, IR decoys 
will have a diminutive effect on the seeker tracking system. 
The radiance of the target is determined within the 
window selected. Values developed by Birk (Ref. 28] for a 
Similar fictitious model with a background of 15°C produce a 
radiance of 2.45 W/sr/m’*. The radiant intensity is determined 
from equation SS Ref. Miz seer cae 
T=(N A) (93) 
where, 
A is the area of the target, mé 
Thus, Upen SubSEM vere y welac- 
di= (2.45) CGea9) 
= 1225 W/sr 


The radiant emittance 1s determined from equation 94 ([Ref. 


12:p. 138}. 
W=t N (94) 
= f- (2.45) 
= 7.7 W/m‘ 
2. IR Range Equation 


From the range equation, the seeker aperture may be 


derived ei. es 2 |Z 
R at OK Ua Ga At Aa WV (95) 
oY [AfAG I” 


where, 


D* is the specific detectivity 


Af is the receiver bandwidth 


WS 


eloecnesarea Of tThe@deeeering material 
CG, is the transmission of the IR optics 
The specific detectivity is a measure of the seeker's 
sensitivity independent of detector material or receiver 
bandwidth [Ref. 27]. If a Mercury-Cadmium-Telluride 
photovoltaic detector (Ref. 27:p. 2.67] is used, then: 
o*% = 10° m-H2"*/ W 
A typical value for the frequency bandwidth is 1 KHz. 
If the desired range, R, is three miles, then aS od 
transmittance of 22% or 0.22. The transmittance of the Opeles 
ls estimated at 90% or 0.9. The area of the detecting 
material will be one square centimeter. [Ref. 12:p. 132] 
Rearranging equation 95 to solve for the aperture 


area, A.,, provides: 


da? 


2 1/2 
A = Ro mw CafAd] 
4 FY Tie Ag W 
Inserting the appropriate values into equation 96 produces: 
é 
Aa= >9.4¢m 
ra 
ees 1 


(96) 


Thus, the aperture diameter is: 


D, = 2.6 inches 
tees moO lanctemasis ea consliaqered reasonable for the 


mrmesile under consideration. 


ea 


C. SEEKER CONSTRUCTION 


The construction of the seeker originates from Wilcox 
(Ref. 30]. The gyro stabilized tracking seeker is mounted on 
a ball support shaft. Two sets of circular ball bearings 
interface between the gyro and the support ball, allowing the 
gyro to spin freely. 

The gyro consists of a ferrous cup, the required optics 
and the detecting material. See Figure 42 for a cross- 
sectional view of the seeker head. 

The ferrous cup is attached to the aft portion of the 
gyro. The outer edge of the cup rotates through the gapger 
four U-shaped coils. Each pair of coils is responsible for 
generating the command torque that will change the gyro's 
Spatial orientation. See Figure 43 for a diagram of the above 
geometry. 

The optics consist of lens that focus the thermal energy 
on to a set of primary and secondary mirrors. The mirrors 
then reflect the target image on to the detecting material. 

The detecting material selected for the window of interest 
is Mercury-Cadmium-Telluride. When this type of material is 


selected, the sensor is referred to as a photovoltaic 


radiation detector. This type of detector works most 
efficiently at 77 K, which implies cryogenic cooling. (Ref. 
2 
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Figure 42. Gyro Cutaway [Ref. 30] 


The detector operates in the following manner. Photons 
emitted from the target enter the detecting material. Upon 


entry of a photon an electric potential is generated. This 


bea, 


potential is measured to determine the amount of thermal 


energy incident upon the detector [Ref. 27:p. 2.58]. 
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Figure 43. Coil Orientation [Ref. 30] 


If four of these detector cells are oriented as shown in 
Figure 44, then four separate voltages can be generated [Ref. 


Beep. 22.87). 
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Figure 44. Detector Cell Orientation 
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D. SEEKER OPERATION 


The missile, in an offensive mode, will be directed to the 
target area by the operator entering the target's range, and 
bearing from the launching platform into the missile interface 
panel. The missile, upon launch, generates internal steering 
commands to position the missile in the approximate vicinity 
of the target. Once the missile is in the target area, the 
seeker is energized and initiates a prescribed track search. 
When the target is located, the detector will attempt to place 
the target image on the optical axis. The optical axis is the 
point where the corners of the four detecting cells meet. See 
Figure 45 for an example of the seeker view. The target image 
is represented by a circle. Cell 3 has the largest portion 
of the image incident upon it; accordingly, the cell generates 
a larger voltage when compared to the voltage output of cell 
10s 

These voltages are processed in the following manner. 
Voltage 3 is added to the inverse of voltage 1. This summed 
voltage is then used to generate a current, i, through coil 
3. The same process is repeated by detector one. In this 
case, the current through coil 3 is increased, producing a 
stronger magnetic field. The increase of the field tends to 
draw the ferrous cup into the gap of coil 3. Likewise, the 


field produced by coil 1 is weakened, allowing the ferrous cup 
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Figure 45. Image Centering Process (Detectors 2 and 4 Not 
Included) 


to withdraw. Thus a couple is produced about the lateral 
axils. The cross product of the angular momentum about the 
optical axis into the couple axis produces a torque vector 
that is oriented downward. This torque causes the seeker head 
to rotate to the right, which consequently centers the target 
image. Once the target image is centered, the voltages are 
in equilibrium and no further torques are required provided 
the image remains centered. [Ref. 30] 

The control system picks off the voltages associated with 
the torque commands and aligns the missile's longitudinal axis 


With the cpticalmaxi1s. 


E. SEEKER EQUATIONS OF MOTION 


This section develops the transfer function relating the 
seeker generated torque to the target error angle. Target 
error angle, ¢~, may be defined as the angle sweep out from the 
optical axis of the seeker to the target's position [Ref. 30]. 

Figure 46 gives the geometry used for developing the 
equations of motion. Wilcox [Ref. 30] gives a complete 
development of the equations of motion. For the purposes of 
this study, the target motion will be restricted to the 
lateral plane of the missile's field of view. Target movement 
in the vertical plane would have to invoke the same response 
as the target movement in the lateral plane because the 


autopilots are identical for these axes. 
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Figure 46. Gyro Coordinates 


For the lateral plane the equation of motion becomes: 
l= J+ a0.¢ (97) 
where, 
W is the nutation angular velocity 
L is the torque required to precess the gyro 
through the error angle. 

Examination of Figure 46 reveals for the optical axis or 
ey 1,, tO precess in the i,, i, piane, a torque must be 
Created about i1,. The magnitude of this torque is given by: 

L=Aw, g (98) 
where, 


A is the moment of inertia about the spin axis 
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We is the angular velocity about the spin axis 

For the steady state case, as the error angle is reduced 
to zero, so should its time rate of change. From equation 98, 
this implies that as the error angle approaches zero, so does 
the torque required. Thus, the torque may be expressed as 
shown in equation 99. 

= -ag+bo(g.-P) (99) 
where, 
a,b are gain constants 
(OoepaO) eS EMG Beis reer epae 

Equation 99 shows clearly that the torque must be filtered 
or damped in some manner to reduce it to zero as the error 
angle approaches zero. 

A block diagram may be constructed as shown in Figure 47 
from equations 6 and 8. Applying the Mason Gain Rule to the 
block diagram, a transfer function for torque/error angle can 
be derived. This transfer function is displayed by equation 
HHO} 6) 


2 
— = _ She a (100) 
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where, 


M = WwW? 


moment of inertia about the 1, 
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Figure 47. 
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Block Diagram of Seeker Equations of Motion 
If the gyro has the properties listed in Table 25, then 
substituting these values into equation 100 yields: 
Boe SSS Se. (101) 
GB 574 524500524 I.xlo’s + 786750 
When a step function is applied to equation 101 and plotted, 
a graph such as Figure 48 is generated. 
The steady state response is of immediate concern because 
"once the steady state is known to be stable, measures can be 
taken to bring the transients into an acceptable level." [Ref. 


me. p. 38) 


3.7 


Table 25. Gyro Parameters [Ref. 19] 
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Figure 48. Seeker Time Response to Step Input (Actual and 
Approximate) 
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Therefore, an approximation of equation 100 may be produced 
to represent the transfer function's response. rie Final 
Value Theorem (FVT) is used to determine the approximation, 
where: 


| a | 
hm f(t) lin Les) o>) 


If ¢ is a ramp function, then equation 102 becomes: 
(S4+m) 

bi ee ok 

20 S°+aust+mset+ bu ~ Ra 





where, 
Ry = mZ (bu) 

Examination of the denominator reveals a pole exists that 
has negligible effects on the steady state solution. Removing 
this pole and adjusting the numerator, yields an approximation 
as shown in equation 103. Applying the FVT to equation 103 


with ¢ aS a ramp input produces: 
(e 
mala Ra Wr S$ (103) 
7 
DO S + Emw,S * We 


4 
lim Ra®n Re b 
a 


G>o 64 LM, S Wy ee i 





Plotting the results of equation 103 to a unit step input 
yields a graph that is similar to Figure 48. Therefore, the 
seeker iS an acceptable design and may be interfaced with the 
autopilot to determine the entire system's time response to 


a given target disturbance. 
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Vill. AUTOPILOT DESIGN 


A. INTRODUCTION 


The arrangement selected for the missile autopilot is to 
use an accelerometer in conjunction with a rate gyro. The 
accelerometer provides a feedback of accelerations developed 
by the missile and the rate gyro acts as a damper. (Ref. 
19:2 1D a 

The scope of the autopilot design approach will be limited 
to the lateral plane of motion. Development of the vertical 
plane autopilot would be very similar to the lateral design 
With the addition of a gravity bias term. This restriction 
also assumes that the roll produced by the missile maneuvering 
is negligible. Justification for the weak roll response is 
well developed in Reference 19. [Ref. 19:p. 92] 

The missile engagement process is initiated when the 
detecting material senses the image of a target. The error 
angle between the target and optical axis is driven to zero 
by the torques created by the seeker system. These torques 
are picked off by the autopilot. The autopilot compares the 
optical axis angle to that of the longitudinal axis angle as 
measured from a common reference. This difference is 


understood by the autopilot as the lateral acceleration 
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required to make the optical axis and longitudinal axis 
collinear. The autopilot generates the command signal that 
allows the missile to achieve the required lateral 
acceleration. This acceleration adjusts the missile's flight 
path so that an intercept with the target is possible [Ref. 


BO]. 


B. STABILITY DERIVATIVES 


Before proceeding with the autopilot design, the stability 
derivatives of the missile need to be resolved. Figure 49 
shows the geometry referenced in the development of these 
derivatives. 

Garnell [{Ref. 19:p. 100] provides a method for estimating 
the derivatives based on desired lateral acceleration. This 
method is accurate up to 0.2 radians of sideslip; thus, the 
missiles motion will be restricted to 0.2 radians. The 
maximum lateral acceleration is assumed to be 32 ft/s‘ with a 
forward velocity of 894 ft/s; Hence, 

y = 22 (104) 
where, ‘ 
Y; is the sideslip angle 
B 1s the normal force coefficient affiliated with 
Sideslip 


If @ = 0.2 radians, then equation 104 becomes 


4e = 160 
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Figure 49. Variable Definition [Ref. 19] 
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For small angles, 


P= V/s (105) 


v is the component of missile velocity along the 


where, 


Y-axis 


u is the missile velocity along the X-axis 


V= (106) 
where Dy Sel” 


VC oeencm ola romeo secoctencilent aaffiliated with 
a component of velocity along the Y-axis 


then substituting equation 105 into equation 106 yields, 


4 = Ka / A (107) 
(160) /(894) 


eC 
Recalling the mass of the missile, 


m, 1S 14.9 slugs and 


the lateral moment of inertia is 11.1 slugs Geo chen che 


coefficient of the moment about the Z-axis for a given 


velocity component along the Y-axis may be determined. 


nyC =Yy mx (108) 
The static margin, x, is assumed to be 5% 


of the missile 
Hemegth or 0.33 ft (Ref. 19:p. 


88] 
Rearranging equation 108 and solving for n 


ives Hy mx / C 


(0.2)(\4,9 (0.3 3)/ CA.) 
0,0% 


yields: 


1" 


y’ 
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The derivatives that involve canard deflections are 
determined next. If the canard is deflected 0.2 radians and 
the canard moment arm (cma) 1S 2.67 feet, then the force 
coefficient due to canard deflection, Yy s 1s determined from 


equation 109. 


Ny UE eee (109) 
= (160) (0.33)/2.67 
= 20 


The moment coefficient due to a canard deflection is 
determined from equation 110. 
Ny “Ym (Cma)/ C (110) 
= 72 
For a canard design this coefficient is positive. 
Garnell [Ref. 19:p. 100] recommends the coefficient of the 
moment about the yaw axis due to an angular rate about the 


yaw axis, n., be given a value of -3.0. This estimate is 


rf 


based on other missile designs. 


C. DEVELOPMENT OF THE AERODYNAMIC TRANSFER FUNCTIONS (ATF) 


There are two ATF of interest for a lateral autopilot 
design. These are the achieved lateral acceleration for a 
canard deflection input and the yaw rate for a lateral 
acceleration input. Euler's equations of motion provide the 
necessary relationships to produce these ATFs. Selecting the 


equations relevant to the lateral plane yields [Ref. 19:p. 
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Sal. 


fyt ver (111) 
V aU Wyo (192 } 
C okallld OS aaa VS Ne% (113) 


From the above equations, a block diagram may be 


constructed as shown in Figure 50. Applying the Mason Gain 
Rule to Figure 50 yields: 
@ 
ty _ Yxo- Yen, S$-Ulng 4, - Ny Yy ) (114) 


2 
A cq Rey TOS ty nt Un, 





Figure 50. {,/%5 Block Diagram 


The remaining ATF is based on equations 112 and 113 and 


the block diagram is shown in Figure 51. MGR produces: 


ei 4 gay Ss Sy (115) 
t S* = (N+ y s+ ny Usen, 4, 





Figure 51. r/z Block Diagram 
The desired result, however, is r/f,. This ATF may be 


realized by inverting equation 114, then multiplying it by 


equation 115. Therefore, 
v : Are str¥y4gny-Ny Yyv (26) 


ty Urs YsNynS- UC Ney, — Ov Ye ) 


D. LATERAL AUTOPILOT OPERATIONS 


The arrangement shown in Figure 52 incorporates a 


accelerometer and a rate gyro. The accelerometer is placed 
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Figure 52. Auto Pilot Block Diagram [Ref. 19:p. 98] 


c feet in front of the CG position and the rate gyro is placed 
such that its position does not coincide with a vibration 
Meee. {Ref. 19:p. 97] 

A description of how the autopilot generates the achieved 
acceleration is as follows. The seeker produces the signal 
that equates to a desired lateral acceleration. This signal 


is fed to the fin servo which generates the proper canard fin 


deflection. The airframe responds to the deflection in a 
Manner represented by equation 11. The airframe is then 
accelerated in the required lateral direction. The 


acceleration produces a body angular rate about the yaw axis. 


The time derivative of this rotational motion is summed with 
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the linear acceleration component accelerometer and amplified. 
The amplified signal 1S combined with the rate gyro output to 
produce a feedback signal. The feedback signal is deducted 
from the desired acceleration signal, thus reducing the error 
signal fed into the autopilot. This process continues until 


the error signal is zero. [Ref. 19) 


E. APPROXIMATIONS OF AUTOPILOT COMPONENTS 


As with the seeker, the initial concern is to examine the 
steady state condition. 
For the fin servo shown in Figure 53, the variables are 
defined as (Ref. 32]: 
K, the motor jgamnpce 2 
a the inverse of the motor time constant, 194 
K. seeker signal gain 
K, forwandfeed galnyeiome 


Kevieloei1ty steedboactaiga ino 0.G 





Figure 53. Fin Servo Block Diagram 
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The values following the above definitions were drawn from 
Vincent [Ref. 32:p. 89]. The fin servo may be represented by 
the transfer function shown by equation 117. 
= | 2 
°/ bare = Kak Ke / (9? 4 (a+ Km kK E+Ka Xe) (127) 
For the fin servo, the natural frequency, WW), Eeumcds fisom - 
2 
= aie! 
bn). Salle ( 


272400 


becomes 


w 
I 


iDiomrad/ & 


The damping ratio, uw, is determined fron: 
= P We 
= 2838 
Rearranging and solving for iu, 


a aoe Ge 2 oC) 


ee 2 


Thus, the fin servo transfer function becomes: 


Tt 


= 22400K,/(9°+ 2684432400) (119) 


REG 
The FVT is used to approximate rye CE 
= Ge = + (il elt ea 120 
ere Un, 4 ea Yoo, Un) ( ) 
= I@l 


then gale (6 represented as a step function, the FVT applied 


to equation 114, yielding, 
in aieve 


2 
eetve = Os ste QI) 


If, 
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a 
then equation 114 reduces to: 


f = e ye a a é 
¥/r = Kew, / (5 +Z2Mse nae 2+ Mrse ee) 
Applying the FVT to equation 19, using a step input, bears a 


Fy (e?* Kae 


Thus equation 122 produces the same steady state value and 


result of: 


is therefore a reasonable approximation to equation 114. 
Garnell [Ref. 19:p. 71] provides an approximation for 


equation 115: 


r Ee eS |) 
pea Va) (ee 


where 


ei 
= “7 (nye Me yy) 


F. AUTOPILOT STABILITY CONSIDERATIONS 


Garnell [Ref. 19:p. 105] derives coefficients that may be 


inserted into the Routh-Hurwity stability criterion as shew, 


below. 


Go, ay CS Cee een 
see 


qay Tr 


(124) 


The coefficients are defined as: 
alae = Ke U ne yy 
ae LS ails 
(\/Ke) + Kglyrr One! 
2M e/ Ke We, 
2 
ae WA oa 


am © 
yp ~ 
oon 
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Selecting the correct gain variable will allow the 
BucoOpilot to satisfy equation 124. For instance, if 
K,= 0.8 
= 7,2) 
K,= 0012 
= at 
as recommended by Garnell [Ref. 19:p. 102} for initial 
estimates, the coefficients become: 
ad, = 10300 
@d = 2160 
ao = 2)4 
a, = 0.74 
GU = 0.004 


Inserting these values into equation 124 shows the autopilot 


mermeet the stability criterion. 


G. 


APPROXIMATION OF THE AUTOPILOT TRANSFER FUNCTION 


If the approximation for the fin servo and the ATFs are 


inserted into the autopilot block diagram, the Mason Gain rule 


may be applied to yield a transfer function relating the 


required lateral 


accelerations. 


in the block diagram. The 


accelerations 


to the achieved lateral 


Figure 52 displays the inserted approximations 


transfer function that can be 


derived from this figure is given by equation 125. 


lived 


fu/ fug = (125) 
4.Sxjo°/ (644291674 S.4x 10%e* + 5.2x10%s+ 4.5x10%) 
The approximation of equation 125 will be based on its 
dominant roots. Garnell [Ref. 19:p. 103] shows a technique 
that can be used to factor the fourth order denominator into 
a pair of quadratics. This technique is highlighted as 
follows. Equate the coefficients of the denominator in 
equation 125 with: 
one Bacne hs + 2095-d=0 (126) 
then solve the cubic, 
2K= bK* + 2( gamed kK = ea tece see _ (227) 
Inserting a real root from the cubic solution set into 
equation 128 and 129 will solve for e and f, 
C= (4742K -b)'* ate | 
f= (kt =d)\/2 (129) 
Substituting these coefficients into equation 130 and 131 


yields two quadratics: 
6°44 (4- 6)S4 (K-f)=0 (130) 


Ses oars (eo EP ee (eae 


Therefore, the denominator becomes, 

(S°e|0SSN21) ( s*igeze+ Sides 
The faster mode, where W, = 226 rad/s, will add but a small 
contribution to the autopilot response and can therefore be 
discarded. As shown in Figure 54, the dominant mode, where 


w, = 11 rad/s, has a response similar to that of equation 125. 
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Figure 54. AutoPilot Time Response 
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This autopilot system can be reduced to a first order equation 
based on the time response. Examination of Figure 54 shows 
a time constant of 0.15 can be used to a generate the first 


order equation, hence: 
f 7 
J . ce 
‘Psy 
H. TARGET ENGAGEMENT SCENARIO 


Garnell [{Ref. 19:p. 203] develops a mathematical model of 
a missile engaging a target. Figure 55 defines the geometries 
involved. The line segment MT, is the original line of Saaiie 
with the target heading 0, degrees from this reference iam 
is a line parallel to MT,; if the target and missile wana 
constant velocity, then the target would be located at Taya 
this case however, the target has accelerated and is displaced 


from the original LOS by @ degrees. This error is equal to: 


d= ARCTAN[(Z,-Z,VR] (132) 
where 
Z,, 4, the distance of the displacement measiaed 
perpendacular to Mi, 
R is the instantaneous range 
The missile must adjust its flight path so that (2. - 25am 


when R= 0 at point I. This adjustment is accomplished by the 


missile producing a lateral acceleration according to: 


fu~ K¢ (133) 


Therefore, this acceleration causes the LOS to rotate to a 


154 





Figure 55. Target Engagement Geometry [Ref. 19] 
feet20n parallel to the original LOS and reduces (2, - Z_) to 


zero. 

Figure 56 shows the block diagram that relates the 
intercept geometry with the calculated transfer function. The 
transfer function 2/2, may be determined from Figure 57, where 
G is the product of the forward feed open loop and H is the 
feedback gain. The results are shown in equation 134: 


OL Cs is ie a>) FOS + (Ot) (134) 


where, 


t= K/ KR 


13> 





Figure 56. Engagement Block Diagram [Ref. 19] 
K is selectable by the designer and R is driven by the 


scenario. If a scenario is created such that an engagement 


— = = —?: i a a 





Figure 57. Transfer Function 


time is 20 seconds, K = 4 and the missile flight path is 
orthogonal to the target's heading. Values for t may be 
produced as shown in Table 25, these values may be inserted 
Into equation ks2- Figure 58 shows the responds of Z to a 
unit step input by Z,. Clearly as the to intercept time, t., 


is reduced the damping on the control system response 


ISS 


decreases. This reduction in damping allows the missile to 
react rapidly to target perturbations. The extreme 
oscillations noted in the figure appear to have made the 
missile unduly stable, but note that t, will equal zero before 
the effect of these oscillations can be exerted on the 
missile's airframe. [Ref. 19:pp. 186-187] 


Table 26. t value vs. Time to Go 


tg (sec) ner le 7s!) 














Figure 59 shows that for a given step input in positional 


feeplacement, 2, the missile will develop a lateral 
Meee rcration to drive (Z, - 2) to zero early in the 
engagement. For example, when t, = 20 seconds, the target 


undergoes a displacement, this causes the line of sight to 
rotate. In response to this rotation, the missile generates 
a small lateral acceleration; the displacement between the 

target and missile along the Y-axis is driven to zero in eight 


seconds after the perturbation has occurred. The remainder 
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Figure 58. Positional Response of Missile 
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Figure 59. Acceleration Response of Missile 
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of the flight would not require any further maneuvers on the 


missile's part, provided the target remains at a constant 


velocity. 


Therefore, the autopilot will allow the missile to track 


and intercept targets of interest. 
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IX. CONCLUSIONS AND RECOMMENDATIONS 


A. OVERVIEW 


The Naval Postgraduate School was requested by the Naval 
Air Test Center to study the feasibility of developing an 
inexpensive, short-range weapon system for use against hostile 
surface combatants. The weapon conceived from this study is 
a 500-pound class, short-range, high-subsonic, anti-ship 
missile. A solid rocket motor propels the weapon to its 
cruise velocity from relatively low launch altitudes and 
airspeeds. See Figure 60 for the an illustration of the 
missile's launch envelope. The 150-pound shaped charge 
warhead is estimated to be effective in neutralizing targets 
up to a displacement of one kiloton. Figure 61 shows a size 
comparison between the missile designed and a current air-to- 
surface missile in operation. 

A unique feature of the missile is its defensive launch 
capability. In this mode, the launching platform can deliver 
an attack rapidly, forcing the hostile target to divert its 
attention from the launching platform. This diversion will 
afford the launching platform time to remove itself from the 


hostile weapons' engagement envelopes. 
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Figure 60. Launch Envelope _ 
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Figure 61. Size Comparison 


The missile uses an inertial/infrared guidance system 
to locate and track intended targets. The missile is guided 
to the target area by a simple onboard inertial system. Once 
in the target area, terminal homing is provided by an infrared 
seeker. 

Another advantage to this system is that the launching 
platform does not require any modification to support the 
weapon. 

Launching data, if employing the OSM, can be entered from 
existing armament systems. Therefore, platform modification 


@@e cs san be eliminated. 


Ge 


The weapon is also designed to produce a low unit cost. 
A conservative estimate of the price per weapon based on a 


limited production of 2,000 missiles is as follows. [Ref. 33] 


Guidance and Comte ee S$ 40 K 
ALYrfYrame..... ses «chee ec 6 ae cee eee S 159m 
MOOR 0 ieee: w e666 seer oe et eee $ 10K 
Warhead s « «6s wea he eee tee reese S$ 5395 
Auxiliary EQULpmen E vere eee S$ Sor 
Total. ..0c +s. > 6 a eee $ 75 K (FY-88 $) 


B. INITIALIZATION AND LAUNCH MODE DESCRIPTION 


The missile has two possible operating modes, which are 
the Defensive Launch Mode (DLM) and the Offensive Strike Mode 
(OSM). The DLM is useful out to 7 miles from the launching 
platform and with a slightly degraded probability of a hit can 
be used out to the weapon's maximum range of 15 miles. The 
OSM is effective out to the weapon's maximum range. 

A DLM scenario may be that a target has become an 
immediate threat to the delivery platform. The pilot, upon 
target identification, aims the nose of the launching platform 
at the target. The weapon station carrying the missile was 
selected prior to entering the operating area. The pice 
moves the master arm switch to the ARM position on the weapon 


armament panel. This switch selection energizes the missile's 
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battery pack, closes the arming fuse circuit, dispenses the 
cooling agent that will bring the detecting material to 77 K, 
and brings the seeker gyro up to operating speed. This 
sequence of events is estimated to take 2 to 3 seconds. The 
pilot waits the appropriate time interval to launch and then 
depresses the WEAPON RELEASE button. The aircraft is free to 
begin evasive maneuvering upon weapon release. The weapon 
free falls from the platform, at which time the missile's roll 
ama pitch are captured. The motor is ignited 1.5 seconds 
after release. At the end of the boost phase, the seeker 
begins to search for the intended target from the release 
altitude. A barometric altimeter is used for pressure 
reference, which allows the missile to maintain its release 
altitude. The barometric altimeter is assumed to be less 
expensive than a radar altimeter as an altitude reference 
device. 

When the target enters the seeker's field of view, the 
missile is given commands as discussed in Chapter 8. The 
missile begins to home on the target's image. Upon contacting 
the target, the shaped charge detonates, releasing a molten 
stream or jet into the interior of the target. 

In the OSM, the weapon is fed targeting data from an 
existing missile data panel. The three-axis autopilot 
accelerometers are used to determine the missile's position 


relative to the target's position. The required steering 
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commands are generated to guide the missile into the target 
area. 

The missile cruises and searches for the target at 500 
feet altitude in this mode. Comparing the barometric 
altimeter pressure at release to the entered altitude from the 
missile data panel and using an assumed pressure differential 
determines the required altitude adjustment to hold the 
missile at the desired search altitude. Figure 62 illustrates 


the flight profiles of the two launch modes. 


C. SEEKER DESCRIPTION 


The seeker is a gyro stabilized Mercury-Cadmium-Tellurate 
detector system. It has the capability of detecting a 500 m*° 
area at a range of 3 miles. The detecting material is divided 
intOntOUnR Bees - When a target is recognized, torques are 
applied to the seeker gyro so that each cell has an equal 
portion of the target image projected upon it. This type of 


cantered alignment will place the optical axis at the center 


of the target. 


D. RECOMMENDATIONS 


To further advance the design of the missile the following 
recommendations are made. A analysis computer code should be 
used to refine the lift, drag, and stability derivatives 


values calculated and to accurately size the lifting and 
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Figure 62. Flight Profile 


control surfaces. With such a code, several flight profiles 
could be examined. 

An examination of warhead size versus range should be 
conducted. For example, if the maximum range could be reduced 
to 10 miles than it may be possible to incorporate a 300-pound 
warhead. 

The control system will require further development beyond 
the lateral autopilot stage. Future efforts might include 
designing the vertical and roll autopilots, designing the 
search altitude hold system, and development of a simple 
inertial system that will be used to guide the missile into 
the target area. 

The structural aspect of the missile was not addressed in 
this study; hence, a preliminary structural analysis should 
be conducted. 

Provided a microprocessor is used to coordinate the 
missile's actions, the software needed to execute the desired 
commands requires development at the assembler lever. 

A thorough cost analysis, including life cycle costs and 
Maintenance, should be conducted to better establish the 


overall program expenditures. 
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